Accuracy and effectualness of closed-form, frequency-domain waveforms for 

non-spinning black hole binaries 



Thibault Damour/'^ Alessandro Nagar/ and Miquel Trias^ 

^Institut des Hautes Etudes Scientifiques, 91440 Bures-sur-Yvette, France 
^ICRANet, 65122 Pescara, Italy 
Departament de Fisica, Universitat de les Illes Balears, 
Carretera Valldemossa km. 7.5, 07122 Palma de Mallorca, Spain 
(Dated: December 16, 2010) 

The coalescences of binary black hole (BBH) systems, here taken to be non-spinning, are among 
the most promising sources for gravitational wave (GW) ground-based detectors, such as LIGO and 
Virgo. To detect the GW signals emitted by BBHs, and measure the parameters of the source, one 
needs to have in hand a bank of GW templates that are both effectual (for detection), and accurate 
(for measurement). We study the effectualness and the accuracy of the two types of parametrized 
banks of templates that are directly defined in the frequency-domain by means of closed-form ex- 
pressions, namely 'post-Newtonian' (PN) and 'phenomenological' models. In absence of knowledge 
of the (continuous family of) exact waveforms, our study assumes as fiducial, target waveforms 
the ones generated by the most accurate version of the effective-one-body (EOB) formalism, cal- 
ibrated upon a few high-accuracy numerical relativity (NR) waveforms. We find that, for initial 
GW detectors the use, at each point of parameter space, of the best closed-form template (among 
PN and phenomenological models) leads to an effectualness > 97% over the entire mass range and 
> 99% in an important fraction of parameter space; however, when considering advanced detectors, 
both of the closed-form frequency-domain models fail to he effectual enough in significant domains 
of the two-dimensional [total mass and mass ratio] parameter space. Moreover, we find that, both 
for initial and advanced detectors, the two closed-form frequency-domain models fail to satisfy the 
minimal required accuracy standard in a very large domain of the two-dimensional parameter space. 
In addition, a side result of our study is the determination, as a function of the mass ratio, of 
the maximum frequency at which a frequency-domain PN waveform can be 'joined' onto a NR- 
calibrated EOB waveform without undue loss of accuracy. In the case of mass ratios larger than 4:1 
this maximum frequency occurs well before the last stable orbit, leaving probably too many orbital 
cycles to be covered by current NR techniques if one wanted to construct accurate enough hybrid 
PN-NR waveforms. This problem will, however, be probably greatly alleviated, or even solved, by 
using the EOB formalism instead of PN theory. 

PACS numbers: 07.05.Kf, 04.30.-w, 04.25.dg, 04.25.Nx, 



I. INTRODUCTION 

During the last 5 — 10 years very significant progress 
has been made towards the detection of gravitational 
wave (GW) signals from binary black holes (BBH) coales- 
cences. On the experimental side, current intcrfcrometric 
GW detectors have been operating at their design sen- 
sitivity for several years [l|-|7[ and advanced versions of 
these detectors (more sensitive by an order of magnitude 
0) are being planned and are expected to be operational 
in four to five years. For example, a (10 -I- 1O)M0 bi- 
nary [Mq denoting the mass of the Sun] can be detected 
out to distances of ^ 160Mpc by initial detectors and 
^ 2200Mpc by advanced detectors 0. On the theoreti- 
cal side, one has recently obtained a good understanding 
of all the stages in the dynamics of coalescing BBHs and 
their emitted gravitational radiation thanks to a combi- 
nation of theoretical techniques: (i) the early, adiabatic 
inspiral stage of the evolution can be described in terms 
of a post-Newtonian (PN) expansio n up to 3.5PN order 
in phase and 3PN in amplitude [l0l - [l3| : (ii) the late in- 
spiral, the merger and the ring-down can be c omp uted 
by numerically solving Einstein's equations |14| - |27| : and 



(iii) the effective-one-body (EOB) formalism [28l - [3]| has 
shown its capability of incorporating information com- 
ing both from PN theory and numerical relativity (NR) 
results into an accurate description [s^ - lioj of the whole 
BBHs dynamics and GW radiation, from the early inspi- 
ral, right across the last stable orbit, to the "plunge" , the 
merger and the final ringdown. 

From the data analysis point of view, the detection 
of the GW signals emitted by BBHs necessitates the 
prior knowledge of 'model waveforms', i.e. of faithful 
representations of the emitted gravitational waveforms. 
These model waveforms are used to detect, and then 
to measure the properties of, a putative BBH by em- 
ploying them to match-filter ^l| the output of the de- 
tector. As the expected GW signals depend on several 
continuous parameters (notably the masses of the two 
BHs in the non-spinning case), one needs to construct 
very large banks of model waveforms, in order to densely 
sample the multi-dimensional space of possible signals. 
Full NR simulations are computationally much too de- 
manding to provide sufBciently long waveforms, for suf- 
ficiently many parameter values. Faster methods, mak- 
ing use of our analytical knowledge of BBHs, need to 
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be considered for data analysis purposes. Several pos- 
sibilities can be (and are being) pursued: (i) one can 
use PN theory by itself to construct template waveforms 
representing the inspiral part of the evolution; and the 
PN templates can be defined either in the time- domain, 
or thejrequency- domain (see (42j for various possibilities 
and [9| for a recent comparison of these possibilities); 
(ii) one can use EOB theory, calibrated to a small sam- 
ple of high-accuracy NR waveforms, to construct time- 
domain waveforms that cover the whole BBH evolution, 
from early inspiral to ringdown; (iii) a third possibility is 
to construct closed- form, frequ ency - domain "phenomeno- 
logical" waveform models [H, lU |93| ; these models are 
obtained by fitting some piece-wise analytical formulae 
to the Fourier transforms of a small sample of hybrid 
time-domain PN-NR waveforms. 

From a data analyst point of view, an attractive fea- 
ture of both the frequency-domain PN models, and the 
phenomenological ones, is that they are defined by closed- 
form expressions directly in the Fourier domain, convert- 
ing the waveform generation to an extremely fast process, 
which is very convenient to perform extensive searches or 
parameter estimation studies. The EOB waveform mod- 
els are obtained, in the time domain, by integrating ordi- 
nary difi^erential equations (ODE's) that depend on the 
continuous parameters of the considered BBH. Each EOB 
waveform can be obtained quite fast, both in the time do- 
main, and then (using a Fast Fourier Transform (FFT)) 
in the frequency domain (see below for numbers). How- 
ever, the exploration of variations in the continuous pa- 
rameters of the source is much more time-consuming than 
in the case of analytic, closed-form frequency-domain 
models. It is therefore useful to compare the closed-form 
model waveforms to the NR-calibrated EOB ones. If one 
were to find that the closed-form frequency-domain wave- 
forms were so close to the NR-calibratcd EOB ones that 
they would be essentially equivalent for performing the 
data analysis of GW signals, there would be no need to 
construct large banks of EOB templates. 

This raises the issue of having quantitative measures 
of the "closeness" of two model waveforms that directly 
translate into their closeness in the data analysis of GW 
signals. We recall that, in GW data analysis, model 
waveforms can be used in two different ways, either for 
the mere detection of the signals, and/or, after detec- 
tion, for the measurement of properties of the source. In 
both cases one uses match-filtering [4l|, i.e. a convolu- 
tion of the measured strain in the detector with a large 
bank of model waveforms. Ideally, one can and should 
(for consistency) use the same bank of waveforms for de- 
tection and measurement. But, as the accuracy require- 
ments are different for detection and measurement (and 
are stricter for measurement), it might be advantageous, 
from a practical point of view, to have in hand two differ- 
ent banks of template waveforms. The issue of defining, 
in a quantitative manner that is related to data analysis, 
the closeness of two model waveforms has been discussed 
in several works (45l - l49| , and has already been applied to 



comparin g di fferent models in, e.g., Refs. I, m, m, [13- 

[53 |. Ref. j45l| introduced a special nomenclature for re- 
ferring to two particular measures of the closeness of a 
model waveform to a supposedly exact waveform: (i) the 
effectualness £ (related to the maximum signal-to-noise 
ratio (SNR) associated to the use of a bank of model 
waveforms); and (ii) the faithfulness J- (measuring the 
SNR associated to a model waveform having the same 
physical parameters as the exact one) . The effectualness 
£ is a useful measure for detection purposes, and we shall 
use it below. The faithfulness T was defined as a quan- 
titative way of gauging the quality of a model waveform 
for measurement purposes. The recent literature j46l l47j 
has emphasized that, for measurement^ there is a bet- 
ter, data-analysis-relcvant, measure of the closeness of a 
model waveform to a supposedly exact one, namely the 
Wiener (noise-weighted) squared distance between the 
two waveforms. Below, after discussing the motivation 
for the use of this measure, we shall introduce a certain 
related quantity that we shall call the inaccuracy I of a 
model waveform (with respect to a supposedly given ex- 
act waveform). Then, we shall use this inaccuracy X to 
gauge the quality of the closed-form, frequency-domain 
model waveforms. [We shall also see that, in most cases, 
the 'inaccuracy' is essentially proportional to the "un- 
faithfulness" 1 — T}\ 

The aim of this paper is to study, by means of the 
two data-analysis relevant measures of effectualness and 
(in-) accuracy, the two classes of closed- form, frequency 
domain waveform models for non-spinning BBH systems, 
namely PN expansions and phenomenological waveforms. 
There remains, however, to discuss the issue of the choice 
of supposedly exact target waveforms. In absence of 
knowledge of the (continuous family of) exact waveforms, 
our study will assume as fiducial, target waveforms the 
ones generated by the most accurate version [13] of the 
effective one body formalism, calibrated upon a few high- 
accuracy numerical relativity waveforms. We recall that 
Refs. [13, [13 have shown that, after calibration of a few 
EOB "flexibility" parameters onto some nonperturbative 
information extracted from high-accuracy NR results, the 
EOB waveforms exhibited an excellent agreement (within 
numerical errors) with state-of-the-art NR waveforms, 
for several mass ratios (1:1, 2:1 and 3:1). One is there- 
fore justified in considering the two-parameter family of 
NR-calibrated EOB waveforms as our current best ap- 
proximation to the (unknown) continuous family of (non- 
spinning) exact waveforms. Independently of this argu- 
ment, it is anyway (as already said above) interesting to 
discuss, in quantitative terms, the closeness between the 
bank of closed-form waveforms, and the bank of EOB 
ones. Given the small computational time required for 
each EOB simulation, we shall be able to exhaustively 
explore the whole parameter space, and to identify the 
regions where the closed-form models lack in effectual- 
ness or accuracy. In the process, we shall also describe 
how to construct frequency-domain hybrid EOB wave- 
forms by joining (in a smooth manner) an early-inspiral. 
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closed-form frequency-domain PN waveform, to a sub- 
sequent Fourier-transformed BOB waveform. This will 
lead us to discussing the maximum frequency where the 
"junction" can take place without undue loss of accuracy. 

This paper is organized as follows: in Sec.|lT]wc briefly 
review the basic concepts and definitions used in GW 
data analysis and discuss the accuracy standards that 
will be used in the rest of the paper. Section IIIII sum- 
marizes the structure of the closed-form waveform mod- 
els considered in our study, mainly with the purpose of 
establishing a common notation. In Sec. IIVI we define 
the time-domain BOB models we use, and describe the 
computation of their Fourier-transform, as well as a pro- 
cedure for constructing hybrid PNUEOB waveforms in 
the Fourier domain. Then, in Sec. |V] we study what is 
the maximum frequency value where the joining (in the 
frequency domain) between the early PN waveform and 
the subsequent BOB one can be done without undue loss 
of accuracy. The main results of the paper, concerning 
quantitative measures of the closeness between closed- 
form models, and BOB ones, are presented and discussed 
in Sec. lVIl Finally, Sec. |VlT] presents our conclusions and 
Sec. I Villi the pointers for future work. 

Also, let us note that we generally use units such 
that G = c = 1. In addition, when dealing with the 
BOB dynamics, and Fourier transform, we shall often 
use (without warning) the adimensionalized time ^eob — 

^physical / • 



II. DETECTING GRAVITATIONAL WAVE 
SIGNALS INTO THE NOISE 

In this section, we start by providing a brief review 
of the basic concepts and definitions used in GW data 
analysis, mainly to establish our notation. Then we give 
the definition of the two quantitative measures ( effectu- 
alness and inaccuracy) that we shall use to quantify the 
accuracy of approximated waveform models in compari- 
son with some (fiducial) exact waveform. All our deriva- 
tions will be done assuming observations from a single 
GW detector, although they could be easily generalized 
to multiple detectors. 



A. Conventions and notation 

The measured strain in an interferometric GW detec- 
tor, s(t), is a time scries data consisting in the combina- 
tion of the detector noise, n,{t), and the exact GW signal 
that we want to detect, hx{t) (here, and in the follow- 
ing, the subscript x stands for "exact"). By assuming 
additive noise, we have 



3{t) = n{t) + h,{t) 



(1) 



Commonly, the noise is modeled to be zero-mean, Gaus- 
sian and stationary, and is characterized by its (one- 



sided) power spectral density (PSD), S'„(/), which is de- 
fined as 



1 



n*{f)n{f')^^SMnf-f') 



(2) 



Here, the overline denotes the ensemble average, the su- 
perscript star indicates complex conjugation, and the 
tilde denotes the Fourier transform. 



n{f) ^ / n(t)e-2"/*dt . 

^ — OO 



(3) 



Note that wc use here the LIGO Scientific Collaboration 
(LSC) sign convention for the Fourier transform, with 
a factor J dte~'^"^* {...). This sign convention is oppo- 
site to the one used in many early GW papers (such as 
Refs. dH, [s^]) which used the (theoretical-physics) 
convention J dte'^^'^^^^ {...). This change of convention 
introduces several sign changes in key formulas, that we 
shall indicate below. [Let us warn in advance the reader 
that, because of this, there are sign errors in some of the 
recent GW literature (such as [55|).] 

In this paper, we shall consider three different ground- 
based detectors: Initial LIGO, Advanced LIGO and Ad- 
vanced Virgo. Following Ref. [11], the one-sided noise 
PSDs of these detectors is approximated by some sim- 
ple analytic expressions of a dimensionless frequency 
X = f/fo (see Fig. [J). For Initial LIGO, /o = 150 Hz, 
and (using a low frequency cut-off, /low = 40 Hz) 



Sn (xif)) = 9 X IQ-'^e [(4.49^)-56 ^ Q lQ^-i- 
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-0.52 + 0.32a;'' 



(4) 



For Advanced LIGO, /o = 215 Hz, and ( using /low — 10 
Hz) 



Sn {X{f)) = 10 



-49 



-4.14 



5x 



111 



(5) 



For Advanced Virgo, /o = 720 Hz, /low = 10 Hz and 
Sn (.t(/)) = 10"^^ [2.67 X IQ-^x-^-^ 

+0.59 e^'" [-3.2-1.08 ln(x)-0.13(ln xf] ^-4.1 
+0.68e-0-^3(ln.)^^5.34l (g) 

Let us note here that the initial LIGO noise PSD cor- 
responds to the one in the Science Requirements Doc- 
ument, which is very close to (but not the same as) 
the actual noise curve obtained during LIGO's S5 and 
S6 science runs. The noise PSDs that wc use for ad- 
vanced detectors are also slightly different to the ones 
currently considered by the LSC and Virgo Collabora- 
tion (see Fig. 2 of Ref. HI). 

See Fig. [T] for a plot of these three PSDs in the 
form of the dimensionless effective GW noise hn{f) = 
(/5„(/))^/^. We recall that /i„(/) is a useful mea- 
sure of the noise in that, together with the correspond- 
ing dimensionless frequency-domain effective GW sig- 
nal hs{f) = \fh{f)\, it yields the squared SNR 
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FIG. 1: Effective GW noise /i„(/) = ifSn{f)Y^^ of the three 
GW ground-based detectors considered in the paper. 



within the detector strain signal s, the latter quantity is 
the sum of a zero-mean random Gaussian variable, say 
N = {hm\n) (filtered noise) and of the non-random "fil- 
tered signal", pm, defined as, 

Because of our use of a normalized template hm it is easily 
checked that the zero-mean random Gaussian variable 
N = {hm\n) has a variance equal to unity. Therefore, the 
filtered signal, pm, represents the mean signal-to-noise 
ratio (SNR) for the detection of a GW signal within the 
detector output s = hx + when using hm as model 
template. The optimal SNR will be obtained when the 
template used is the exact signal that wc want to 
detect, i.e. 



as an integral over the loga rithmic frequency: p^ = 

jr{df/mAf)/hM)?m 

The assumed Gaussianity of the noise naturally leads 
to the definition of an inner product (Wiener scalar prod- 
uct) between two (real) time series, a{t) and b{t), namely 



(7) 



This Wiener scalar product is real and symmetric. The 
associated norm, say 



(a|a) = 4 Re / d/ 



a*{mf) 

Snif) 



(8) 



is positive definite, and endows the space of (real) signals 
with an Euclidean structure. 

Moreover, the normalization of this norm is such that 
the probability of having the realization of a particular 



noise series n{t) is p{n) 



From this defini- 



tion, the Neyman-Pearson criterion leads to an optimal 
search statistic equivalent to the form {hx\s) (Wiener fil- 
ter), cf. Eq. (A24) of Ref. [S^. Here, we are using a 
hat to denote normalized templates h = 1 so that 

{h\h) = 1. By analogy with the usual Euclidean space, 
the inner product between two normalized signals, say 
varies between —1 and -1-1 and can be thought of as 
defining the cosine of the "angle" a{h,g) between the two 
corresponding "unit vectors" . One often refers to this co- 
sine as being the "overlap" between the two normalized 
signals, say 



0{h,g) = cosa{h,g) = {h\g) . 



(9) 



As we do not know the exact signal h^, but only some 
approximate model of it, say hm (where the subscript m 
stands for 'model'), one uses as practical search statis- 
tic the quantity {hm,\s). In presence of a non-zero signal 



Popt = {K\K) - = ihxih^y/' . (11) 

In general, we can represent a complex gravitational wave 
signal in the Fourier domain in its polar form as 



Hf) = Ah{f)e-'^-^f^ 



(12) 



where the amplitude Ah{f) and the phase iphif) sue real 
quantities. Since we are considering here real signals, we 
have h{f) = h*{—f), which implies that the amplitude 
and the phase are, respectively, even and odd functions 
of the frequency. Note that here we have explicitly in- 
troduced a minus sign in the phase to ensure, in view 
of the (LSC) sign convention Eq. (jS]) in the definition of 
the Fourier transform, that tpi^) (with lu = Stt/) is, in 
the stationary phase approximation (SPA), the Legen- 
dre transform of the time-domain phase (modulo —j)- 
ip{uj) = ujt — (t>c,w{'t) — J, with d4)c,vi{t) / dt = w . This 
has the advantage that the so-defined frequency-domain 
phase ■)/'(/) then coincides (possibly modulo— -j) with the 
SPA phase tl^if) used in many other papers (such as 
[4^ [50|). The reader should be wary of possible errors 
in the literature due to the change in the sign convention 
for the Fourier transform, as many early papers used the 
theoretical-physics convention. 

Note, for future use, that, in terms of real amplitudes 
and phases, the Wiener products {h\h) and {h\g) read 



{h\h) 
{h\9) 



= 4 



d/ 



d/ 



Snif) 

A^if)Agif) cos iMf)~Mf)) 

Snif) 



(13) 



Before proceeding to the definition of the accuracy 
standards of a particular waveform model, let us spell 
out the notation we shall use for the different waveforms 
that will enter our analysis: 
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hx Exact waveform 

hjn Model waveform having the same physical pa- 
rameters as hr,. 
hfh Model waveform having parameters that max- 
imize the overlap with hx 
hi, Model waveform of a certain discrete template 
bank having a maximum overlap with hx- 
Finally, let us mention that we shall be neglecting here 
the eventual calibration errors from the detectors. Tak- 
ing them into account would imply a strengthening of the 
accuracy standards discussed below (see pTf). 



B. EfFectualness of model waveforms 

Any modeled waveform h„i represents an approxima- 
tion to the real gravitational wave pattern hx that is emit- 
ted in nature, so that the practical search statistic (ft.,n|s) 
will necessarily be less efficient in detecting a signal than 
the optimal Wiener filter {hx\s). In GW searches, one 
says that hm is an effectual model of the exact signal if it 
allows for a successfuU detection of the presence of a GW 
signal in the detector output, even if the values of param- 
eters entering h^ are not close to the real values of the 
parameters entering the exact signal hx ■ [One sometimes 
even allow the parameters of hm to take physically un- 
realizable values.] Following [4^, we shall quantify the 
effectualness, say £, of a continuous (multi-parameter) 
bank of model waveforms by considering the maximum 
value of the ratio between the SNR, pm, obtained by 
using hm as filter, and the optimal SNR, popt = Px, ob- 
tained by using the exact Wiener filter {hx\s), i.e. 

effectualness = £ = {hm\hx) = , (14) 

Popt 

where we recall that hm stands for the normalized wave- 
form in the continuous bank of models that maximizes 
the overlap with the exact (normalized) signal. In other 
words, if we represent the set of parameters that charac- 
terize our model by the vector A; hm is such that 



{hx\hm) ^ ITijpi {hx\h^^^^) 



(15) 



In geometrical language, this means that the cosine of the 
angle (in Wiener's signal space) between the unit vector 
hx and any unit vector in the continuous bank h^^^-^ is 

maximized^ by the unit vector hfn- This maximum co- 
sine, corresponding to a minimum angle, will necessarily 



range between —1 and -1-1, and in practice will be close 
to, but strictly smaller than, -1-1. 

In this paper, we will restrict the source parameters to 
take values with real physical interpretation. In partic- 
ular, the symmetric mass ratio, v = mim2/{jni + '712)^, 
will be constrained to be < 0.25. [Otherwise, we would 
be considering systems with imaginary individual masses, 
mi and m2 !] 

It is also convenient to work with the ineffectualness, 
i.e. the quantity complementary (with respect to 1) of 
the effectualness, say £ = ! — £. In other words, the 
ineffectualness of a continuous bank of models hm in de- 
tecting a certain exact signal hx is given by 



£xvri 1 {hx\hm. 



(16) 



In choosing an upper bound on the value of the inef- 
fectualness one must, however, remember that, in prac- 
tice, one approximates the continuous bank of model 
templates by some finite, discrete template bank. The 
closeness of a discrete bank to its continuous version is 
measured by the so-called minimal match (MM) (soj . i.e. 
the minimum overlap between any model hm hi the con- 
tinuous bank, and any model /if, in the discrete bank. 
One can replace the minimal match by its complemen- 
tary, i.e. the maximal mismatch smm = 1 — MM. It is 
easily seen (by viewing the complementaries of overlaps, 
1 — cos a ~ as squared distances between nearby 

points on the Wiener sphere of unit signals) that, when 
setting thresholds on the ineffectualness of the discretized 
version of a template bank (say, b), one should replace 
the continuous ineffectualness £xm by the "effective" in- 
effectualness 47 



£xb — £xm + ^MM 



(17) 



Since in any GW observation, the distance to the source 
is inversely proportional to the observed SNR and the 
number of events is proportional to the observable vol- 
ume, we can express the fractional loss of potential events 
in terms of the effective ineffectualness of a given discrete 
template bank. 



# events b ^ n3 

■ = M — txb) 



# events 



total 



(18) 



A threshold on the allowed fraction of missed events will 
set up the effectualness condition for any approximated 
template bank. Typically, one works with template 
banks with a maximal mismatch ^, emm, of 0.005 — 0.030 
[60l - [63j and allow for a fraction of missed events smaller 
than 10% {£xb < 0.0345), which making use of Eq. p7|. 
translates to the following effectualness condition, 



£xrn = {hx\hm) > 0.9705 - 0.9955 



(19) 



Rigorously speaking we should be talking about a supremum 
rather than a maximum, as there is no mathematical guarantee, 
for a continuous bank, that the supremum is reached for some 
model. Anyway, in practice one always uses discrete banks where 
the maximum will be attained by some model. 



^ The choice of this value depends on the computational cost to 
generate the resulting total number of templates and the maxi- 
mum false alarm rate allowed for the search. 



6 



In view of the fact that it is much more difficult, and ex- 
pensive, to improve the sensitivity of the detectors than 
the quahty of the template banks, one should aim at 
imposing on £xb a bound significantly stricter than that 
corresponding to a loss of 10% of the hard-won potential 
events. Then, if we assume the reasonable condition that 
the two contributions on the r.h.s. of Eq. (jl7|) are com- 
parable, i.e. eMM ~ £xm, we see that we should really 
aim at imposing bounds on the effectualness which are 
on the stricter side of the condition above. 

Let us end this section by a brief discussion of the 
procedure we used to compute the effectualness, i.e. to 
maximize the overlap {hm\hx) over all the parameters 
entering the model waveform hm- There are two types 
of parameters we should maximize upon: (i) the physi- 
cal parameters of the source (which, in the present case, 
are the two masses mi,m2, or equivalently the total 
mass M = mi + 7712 and the symmetric mass ratio 
V = mim2/M^); and (ii) some kinematical parameters. 
In the present case, the relevant kinematical parameters 
of hm arc some time origin i^^, and some initial phase 

Among these parameters, the only one upon which 
we can maximize analytically is the arbitrary phase 
Note first that, in the following, we shall work with a 
given target waveform (including a quasi-random ini- 
tial phase), so that our problem here is only to maximize 
over the phase of the template waveform hm. [See ap- 
pendix B of [4^ for the case where one varies the ini- 
tial phases of both the template and the target.] To do 
this we start by noting a simplifying feature of the fre- 
quency domain, with respect to the time domain. In the 
frequency domain, shifting the phase of any waveform 
h{f) = Ah{f)e~'^^'^''^^ simply corresponds to adding an 
arbitrary constant, say 0o, to the phase iphif)- A look 
at the explicit form Eq. of the Wiener scalar prod- 
uct shows that two waveforms whose frequency-domain 
phases tph{f), ipgif) only differ by n/2 are orthogonal 
to each other. [This property does not hold in gen- 
eral when shifting time-domain phases by 7r/2; except 
when: (i) one considers just one harmonic, and (ii) 
the signal is such that the stationary phase approxi- 
mation is accurate, see Sec. IIII Al ] In addition, two 
waveforms having the same frequency-domain amplitude 
Ah{f) have the same Wiener squared length {h\h). This 
shows that the two model waveforms {0°^^ ,0°='^^^} 
define an equal-length, orthogonal basis of the two-plane 
spanned by all the phase-shifted models /i*^". More- 
over, a generic phase-shifted model can be written as 
/i-^o ^ cos(l)oh'l'°=° + sm(j)oh'^°=''/^. These properties al- 
low one to compute the analytical dependence upon 0o 
of the scalar product of h'^" with any given waveform g: 

{g\h'^°) = cosM9\h")+smM9\h''^^) (20) 

With this, one can maximize {g\h'^'') analytically over (po, 
with the result 



On the other hand, the maximization over the other 
parameters must be done numerically. There is, however, 
a simplification when dealing with the maximization with 
respect to the time origin t^- Indeed, the dependence on 
the time origin is known analytically, so that there is, in 
particular, no need to recompute the full waveform every 
time we need to shift a time origin. Indeed, we have for 
the phase difference between any two waveforms 

Mf, to. h) - Mf^ to,a)^ Mf, 0) - i^gif, 0) + ■ 

(22) 

where Ato = to, ?i ~ ^0, g- Then, it is just a matter of 
numerically maximizing the integral of {h\g) given by 
Eq. Uni) over Ato. 



C. Accuracy standards for a model waveform 

Gravitational wave observations will open a new obser- 
vational window on the Universe. In order to do actual 
astrophysics with GW, we need not only to detect sig- 
nals, but also to extract reliable physical information out 
of them. In other words, we need to have at hand wave- 
form models that are close enough to the exact signals to 
allow for an accurate enough measurement of the physical 
properties of the GW source. [As mentioned above, here 
we could be talking either about the same template bank 
used for detection, or about a complementary template 
bank used for analyzing in detail signals having passed 
a first detection threshold.] It has been recently empha- 
sized (iii - fist that a simple accuracy standard ensuring 
that a model waveform hm can be used for extracting 
most of the scientific content of the GW signal h^ (buried 
within the detector noisy output .s ~ h^ + n) is that the 
Wiener squared distance between hm and h^ be smaller 
than one: 



\hm-hx\' <\n\^ = l 



(23) 



where the over line denotes the ensemble average. Here 
it is understood that h^ and hm share the same physical 
parameters. [But, one is allowed to vary purely kinemat- 
ical parameters such as an initial time a-nd an initial 
phase (pQ.] 

There are several ways of thinking about the mean- 
ing of the inequality Eq. (^5)) . Ref. [l^l motivates it 
by considering the accuracy with which one can mea- 
sure the parameter ^ gauging the affine distance along 
a straight line connecting h^ to hm in Euclidean Wiener 
space: h{^) = {l — ^)hx+S,hm- Assuming Gaussian noise, 
the variance of the measurement of ^ is found to be given 
by 



CT? = \h,n — hr, 



max(ff|/7*o) = hg\hy ^ ^g\h^n f 



1/2 



(21) 



(24) 



Therefore, the inequality (j23p implies that the error in ^ 
is larger than its variation when it interpolates between 
the two models, > 1. 
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Another way of motivating Eq. ([23]) is to think in terms 
of the Wiener geometry of signal space. If the detec- 
tor were noiseless, a measurement of the detector output 
s would give us the exact location, in signal space, of 
the incoming signal h^, namely = s. However, the 
presence of noise in the detector is introducing a certain 
uncertainty in our knowledge of h^- More precisely, un- 
der our assumption of Gaussian noise, the probability 
distribution, in signal space, of the measured signal s is 
p{s) (X e^al*^''^!'. Therefore, at, say, the 68% proba- 
bility level, our knowledge is only that the exact signal 
hx lies somewhere within a "unit ball" in Wiener space, 
centered around the observed signal s: \hx — s\'^ = l. The 
maximum knowledge we can extract from the GW data 
about the parameters A of the source is therefore that 
\hx{X) — s\ = \n\, with |np = 1. However, we do not know 
the exact form of the functional dependence of the sig- 
nal upon the source parameters , i.e. the function hx{X). 
We only know an approximation of this functional de- 
pendence, namely the one given by the waveform model 
/im(A). By using the triangular inequality in (Euclidean) 
Wiener space, we can write the inequality 

\h,n{X) -S\< \h,n{X) - hx{X)\ + \hx{X) - S\ 

= IK-^iX) ~ hxiX)\ + \n\ . (25) 



This inequality shows that the condition Eq. ([23| is such 
that, in the worst case"^, it degrades our knowledge on 
the values of the source parameters by being equivalent 
to doubling the noise level of the detector. Let us recall 
that our only knowledge about the source parameters is 
derived from the fact that the measured signal s satis- 
fies \hx{X) — s\ ~ \n\, so that it is the noise level in the 
detector which determines the statistical errors in pa- 
rameter estimation. Therefore, a doubling (in the worst 
case) of the noise will add systematic biases to measured 
parameters of the same order as the statistical errors. 
This is clearly the minimal requirement we should put 
on the difference between km and h^. However, as in the 
discussion above of the threshold for effectualness, we 
should remember that it is much more difficult, and ex- 
pensive, to improve the sensitivity of the detectors than 
the quality of the template banks. Therefore, we should 
not allow the lack of accuracy of waveforms to lead to any 
significant (effective) increase of the detector noise level. 
In other words, we should strengthen the condition ([23)) 
into a condition of the form 



\h^{X) ^ hx{X)\^ < e^\n\ 



2 = e2 



(26) 



where e < 1 must be chosen so as to limit the effective 
increase of the noise level on the r.h.s. of the triangular 



^ As we are discussing the case where we have only one observed 
datum s, one cannot use the fact that "on average" s is centered 
around hx to try to get a better general inequality on the distance 
between h™ and s. 



inequality Eq. ([25]) to an acceptable level. Note indeed 
that the r.h.s. of Eq. ()25p now reads < (1 -I- £)\n\. One 
should take at least e ^ 1/2, and probably e ^ 1/3 or 
more^. Only with such small values of e can we really 



consider /i„ 
hx- 



as being effectively indistinguishable from 



D. 



Definition of an inaccuracy functional 



It is useful to rewrite the accuracy condition Eq. (^5)) in 
a different form, which highlights its connection with the 
faithfulness functional, introduced in (45l |. The l.h.s. 
of Eq. ([26]) depends on the absolute magnitudes of the 
signals (and therefore on the distance to the source). It 
also depends on the normalization of the Wiener metric, 
i.e. on the overall normalization of the PSD of the noise. 
Let us instead define the following functional of h„i and 
hx- 
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\h^iX) ^ hxiX)\^ 



(27) 



We shall call hx] the inaccuracy of the model wave- 

form hm. with respect to hx^- The functional T[hm.;hx] 
remains invariant both when one multiplies hm and hx 
by a common factor, and when one multiplies the noise 
density S'„(/) by an arbitrary factor. In other words, it 
depends only on the shape of the noise PSD, and on the 
relative shapes of the two waveforms hm and hx- As ex- 
plicitly indicated, both waveforms refer to the same val- 
ues of the physical parameters of the source. However, 
we can still maximize over the kinematical parameters 
of the waveform model km, namely some arbitrary time 
origin t^, as well as some initial phase t/)^. 
Using the identity 



|a-5p _ \ct\ 



\b\' 



\b\ 



2 (1 - cos a) , (28) 
\b\ 



valid in any Euclidean space (with a denoting the angle 
between the two vectors), we can write more explicitly 



^ Note that, according to Ref. [64j . the consideration of calibration 
errors in GW detectors also lead to suggesting a strengthening 
of the original condition Eq. I I26I I by a factor tcalib smaller than 
one. But such a 'calibration factor' ecalib < 1 should be taken 
on top of the one of Eq. I I26I1 . which has a different motivation. 
In other words, one should use a "safety factor" etotal = ^ ' ^calib- 

^ Note that this inaccuracy functional is not a symmetric function 
of its two arguments. One could also consider defining the sym- 
metric functional obtained by replacing in Eq. I I27I1 the denomi- 
nator \hx(\)\^ by \hm(^)\\hx(\)\. However, this would lead to an 
upper limit on this symmetric functional which differs from the 
inverse squared SNR by a factor involving the ratio l/iml/l/ia:!. 
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the inaccuracy functional in the foUowing expanded form 



1/2 



{Kn\h„i) 



+2 



1 - {K\K-,) 

(29) 



With this notation, the accuracy condition Eq. (|26p can 
now written in the fohowing form 

^2 



T[h,n\ hj] < e 



{n\n) 



(30) 



where = = denotes the squared SNR. 

This reformulation of the accuracy condition has several 
useful features: 

• All the dependence on the absolute magnitudes of 
both the signals and the noise is rejected to the 
r.h.s., and in a very simple form involving the 
(squared) SNR. This allows one to compute the 
'inaccuracy' on the l.h.s. from the shapes of the 
signals and the shape of the noise PSD indepen- 
dently from the various accuracy thresholds set by 
the r.h.s. 



• The inaccuracy I[/im; ^x], Eq. (|29p . is seen to con- 
sist of the sum of two positive-definite contribu- 
tions: (i) the first one depends only on the ra- 
tio of the Wiener squares of the two signals (and 
vanishes when |/i,„p = while (ii) the sec- 
ond one is proportional (and nearly equal) to twice 
the 'unfaithfulness', ITmx = 2(1 — J^mx), where 
J^mx = {hx\hm) [Hi- Note that, using the explicit 
frequency-domain expression of the inner product 
in terms of the amplitudes and phases of the wave- 
forms, Eq. (jl3p . the first contribution to the inaccu- 
racy only depends on the frequency-domain ampli- 
tudes, Am{f) and A^if), and vanishes when they 
coincide; whereas the second contribution, propor- 
tional to 2Tr,ix, vanishes when both the amplitudes, 
and the phases tjjrnif) and ipxif) coincide. 

• Note in passing that if, instead of evaluating 
T[hjn] hx] for the same values of the physical param- 
eters, one were considering the analog of the 'efFec- 
tualness' f, i.e. if one were minimizing X[hm]hx] 
over the physical parameters of the model wave- 
form, the minimization over the effective distance of 
hm would require the condition = \hx\'^ (be- 
cause of the positive-definite character of the first 
contribution (i)), so that T[hfa \ hx\ would reduce to 
the 'inefi'ectualness' 2(1 — £) = 2(1 — {hx\hfn))- 

• In our comparisons below, our experience has been 
that the inaccuracy contribution from the ampli- 
tude differences was typically several orders of mag- 
nitude smaller than the one from the phase errors. 
This means that the accuracy condition Eq. ([26| is 
practically equivalent to the following faithfulness 
condition 



1-^=1- {hx\hm) < 



2p2 



(31) 



If (as in [43,14^) we expand the inaccuracy to lead- 
ing order in the amplitude and phase differences, 

SA{f) = AraU) - Axif), S^'if) = V™(/) - Mf), 

it becomes 



' SA 

X[hn,]K] = (( ^ 



) + {5^) 



(32) 



where the angular brackets denote an average with 
respect to the normalized measure entering the 
squared SNR, i.e. dwif) = 4iVd/A,(/) V5„(/), 
with the normalization factor N = 1/ p^p^ chosen 

so that dw{f) = 1. This expanded expression is 
too coarse to be numerically useful, but it is useful 
in indicating that the accuracy condition Eq. (|26p 
essentially means that 



6A 



(33) 



i.e. that the noise-weighted quadrature sum of the 
fractional amplitude discrepancy and of the phase 
discrepancy (between hj^if) and h^^f)) must be 
smaller than e/p. [As we shall further discuss be- 
low, this motivates one to considering pcff = p/e 
as an effective SNR.] If, for instance, we consider a 
SNR of 10 and take e = 1/2 (corresponding to an 
effective SNR pcs = 20), wc see that the accuracy 
condition is restricting the noise-weighted quadra- 
ture sum of the fractional amplitude discrepancy 
and of the phase discrepancy to be smaller than 
0.05 (i.e. 5 % or 0.05 radians). Such a constraint 
is more difficult to satisfy for the phase (for which 
any error in its dynamical description is accumu- 
lated over many GW cycles). 

In summary, we shall require of any model waveform 
that it satisfies the following conditions: (i) for detection 
the effectualness condition Eq. p^ : and (ii) for mea- 
surement the accuracy condition Eq. pop . where the in- 
accuracy functional is defined by Eq. ([27]), or, in more 
explicit form, by Eq. ((29|). From its r.h.s. we see that 
the higher the observed SNR, p, is, the more restrictive 
the accuracy condition becomes. In addition, one must 
choose some e < 1 to guarantee that the inaccuracy of 
the model waveform only increases the detector noise by 
a factor 1 + e remaining close to 1. 



III. CLOSED-FORM, FREQUENCY-DOMAIN 
WAVEFORM MODELS CONSIDERED IN OUR 
STUDY 

In this paper we consider GW observations from a sin- 
gle interferometer, assuming the detected signals to be 
short enough so that the orientation of the detector (due 
to Earth's motion) can be considered constant during the 
time it is being observed. [For simplicity, we shall also 
not explicitly write down the (Doppler) effect on the GW 
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signal coming from the relative motion (and cosmologi- 
cal redsliift) between the detector and the source.] These 
approximations arc commonly used for first and second 
generations of ground-based detectors, and will allow us 
to reduce all the dependence on the sky location, lumi- 
nosity distance, polarization and inclination angles to an 
effective distance parameter, D^s (sec Eq. (|36p below), 
and to an additional phase shift ipo (sec Eq. (|37p be- 
low). Moreover, we shall only consider the dominant 
(2, 2) mode from the spin-weighted spherical harmonics 
decomposition of a signal from a non-spinning coalescing 
binary system. 

In a suitable frame adapted to the quasi-circular orbit 
of the considered compact binary system, the two orthog- 
onal polarizations of the emitted time-domain quadrupo- 
lar GW signal read 



h+it) = 



1 



aGw{t) cos (^Gw(i) : 



hx{t) = [cos l] (<) sin (/)gw {t) 



(34) 



The relative weight between the polarizations depends 
on the inclination angle, t, i.e. the angle between the 
angular momentum of the system and the propagation 
direction. The time series acw{t) and (pewit) define the 
(time-domain) amplitude and phase of the quadrupolar 
waveform. [The time-domain phase (pewit) is convention- 
ally defined so that it increases with time] 

The corresponding GW signal seen in the detector will 
be given by a single time series strain, hit), resulting 
from the projection of the two GW polarizations with 
the beam antenna patterns, F+ and Fx , which depend 
on the sky location of the source (i.e. the relative orien- 
tation of the line of sight and the normal direction of the 
detector's plane) and the polarization angle [g^. Using 
basic trigonometric relations, the detected strain can be 
written as a single cosine function, with amplitude a{t) 
and phase (pit), as follows: 

hit) = F+h+ +Fxhy = Cacwit) cos(^c.w(0 + fo) 
= ait) cos (Pit) , (35) 

where we defined 



C = -y'(l + cos2t)2F2 +4cos2iF2 (36) 



^0 



arctan 



—2Fx cos L 
F+(l + cos2 l) 



(37) 



If we were considering observations from multiple detec- 
tors, the amplitude and phase differences between the 
various detected strains would involve different C and ipo 
values through the antenna patterns of each detector, 
and thereby we could infer the sky location and polar- 
ization angle of the source. Also if we were observing a 
GW signal for a long time period, the time dependence 
of these three parameters could be used to estimate the 
position of the source. In this paper, we are considering 
short-lived signals detected by a single interferometer and 



therefore C and cpo will be numerical constants and ab- 
sorbed, respectively, into an effective distance and the 
initial phase of the signal. 

In the Fourier domain, the measured GW signal can 
be represented generically in the polar form Eq. (|12p . 
In what follows, we shall describe the different waveform 
models considered in our study in terms of the frequency- 
domain amplitude and phase, Aif) and ipif)- 



A. Frequency-domain Post-Newtonian model: 
PN(f) 

In the limit where the two coalescing compact objects 
are still far away from the merger, the post-Newtonian 
approximation to Einstein's theory provides an accurate 
representation of the dynamics and radiation of the sys- 
tem. For non-spinning black holes binaries, the time- 
domain phase, (pp^it), is known up 3.5PN order and the 
amplitude, ap^it), up to 3PN pi]| - [l3j . 

A closed-form expression for the PN waveform in the 
Fourier domain can be obtained by making use of the 
Stationary Phase Approximation (SPA). Using this ap- 
proximation is not an extra assumption beyond the usual 
ones underlying the PN method. Indeed, the derivation 
of the time-domain PN phasing already assumes that 
the inspiral is adiabatic, which is the basic assumption 
needed for using the SPA. More precisely the SPA is 
valid when the amplitude and frequency of the signal 
evolve on time scales much longer than the orbital pe- 
riod: Ea = ait) / iait)^it)) < 1 and e/ = 0(t)/^2(t) < i. 
For a coalescing binary system the time scale of evolu- 
tion of the GW amplitude and frequency is the radiation- 
reaction time scale so that ea ~ e/ ~ [v being the PN 
expansion parameter v = (tt/M)^/'^] stays much smaller 
than one essentially up to the end of the inspiral. See 
[sol . l66j for a more detailed discussion of the applica- 
bility of the SPA, and for the calculation of the first 
correction to the leading-order SPA. Here we shall use 
the leading-order SPA which has the following expres- 
sion in terms of the time-domain GW amplitude a(t), the 
GW phase (pit), and the derivative of the GW frequency 
Fit) = 0(t)/(27r): 



hif) 



<tf) 



,-i[27r/t/-f-0(t^)] 



(38) 



Here tf is the saddle point defined as the instant when 
the GW frequency Fit) is equal to the Fourier-domain 
argument /. It is obtained by solving the equation (pit) = 
27r/ in t. 

By inserting in the general SPA expression (|38|) the 
time-domain PN expressions one can obtain a closed- 
form expression for the PN-approximated frequency- 
domain inspiral waveform. As discussed in [4^ there are 
actually several different ways of using the PN informa- 
tion to define hif). Here, we shall consider the most 
straightforward way of doing so (called 'TaylorF2' in 
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[42|). consisting in expanding both the frequency-domain 
amphtude and phase in powers of the PN expansion pa- 
rameter v{f) = {nfMY^^. We recall that M = nii +m2 
denotes the total mass of the system (mi and m2 being 
the individual masses), and v = mim2 / M'^ the symmet- 
ric mass ratio. In the following, we shall denote v{f) 
simply as u, but one must remember that it is a func- 
tion of the Fourier variable /. For brevity, we shall refer 
in the following to this PN-expanded frequency-domain 
waveform as 'PN(f)'. We can write it as 

KM) = ^PN(/)e-'*-(-^'0(z;.so - v) , (39) 

where we introduced a frequency-domain cut-off at the 
Schwarzschild-like Last Stable Orbit (LSO) frequency 
?;(/lso) = 6~^/^ (^(2;) denoting Heaviside's step func- 
tion) . 

The angular argument i/jpnif) hi Eq. (j39p is given by 
the PN series 

where tc and 0c represent some arbitrary time and phase 
shifts. [We absorbed the SPA phase shift —j into f/fc-] 
The explicit expressions of the PN expansion coefficients 



at (up to 3.5PN) can be identified from Eq. (3.18) of ^. 
[Notice that Ref. Q uses the opposite sign convention 
for the Fourier transform, which translates into a global 
sign in the frequency-domain phase ipif), usually written 
explicitly in the definition of h{f) so that ■(/'(/) has the 
same expression with both sign conventions; compare, for 
instance, Eq. ([12]) with Eq. (3.17) of 

In many early GW papers, the GW amplitude was 
taken only at leading order in the PN expansion, i.e. at 
the 'Newtonian' (frequency-domain) approximation: 

Mf)^-Jf/-^v-y^. (41) 

Here we have kept an explicit negative sign coming 
from the second derivative of a complex exponential func- 
tion, and we have defined the effective distance Dcs = 
Dj^/C which absorbs all the sky position, polarization 
angle and orientation dependencies, see Eq. (|36|) above. 
The 3PN correction to the (quadrupolar) frequency- 
domain amplitude is then obtained by combining the 
3PN-accurate expression of h'^^{t) with the 3PN ex- 
pansion of F (given, e.g. , in Q). Re-expanding the result 
as a single PN series leads to 



^pn(/) - AM) 



1 



451;^ 



168 
34 473 0791^3 

6 386 688 ' 



323 

224 



105 271J/2 1975 055:/ 27 312 085 



/857r 



{Aiy - 1) - 24w 



24192 338 688 8128 512 7 V 64 

3 248 849 057^2 545 384 828 789i/ 2057r2t/ 177 520 268 561 428i7r\ . 

1 \ 

105 / 



178 827 264 



5 007163 392 
I 



48 



8 583 708 672 



(42) 



Note that, contrary to the convention generally assumed 
in the rest of the paper, this '3PN amplitude' Ap^{f) 
contains some (small) imaginary contributions (as h^^ (t) 
in [121) ■ When, in the following, we shall talk about the 
amplitude of the PN(f ) waveform we shall mean the mod- 
ulus |ylpN(/)|. Correlatively, when we shall talk about 
the phase of the PN(f) waveform wc shall mean (minus) 
the argument of the complex number /ipn(/), defined in 
Eq. ([39]) (i.e. '0pn(/) minus the argument of Apx(/)). 

Contrary to the time-domain amplitude |ft.(i)|, which 
is dimensionless, the amplitude of a Fourier-transformed 
GW signal A{f) = \h{f)\ has the dimension of a time. 
By contrast, the effective GW signal amplitude hs{f) = 
\f\A{f) = \fh{f)\ is dimensionless. As dimensionless 
quantities, the effective amplitudes of the signal, hs{f), 
and of the noise, hn{f), are not affected when one rescales 
the frequency, by means of the total mass M, into the 
dimensionless variable fM. As mentioned above, the 
SNR2 is given by the logarithmic frequency integral: 
= iridf/f){hs{f)/hn{f))^. The effective GW am- 
plitude of the PN(f) waveform is represented (for the 



equal- mass case ;/ = 1/4) in Fig. O This figure also plots 
the effective GW amplitudes of the other waveforms we 
shall use in the present study: the EOB one, and the 
two phenomenological ones that we shall introduce next. 
Fig. [5] represents also the difference between the phase of 
the PN(f) signal, and that of the corresponding EOB(f) 
signal. This figure indicates that both the amplitude and 
the phase of the PN(f ) waveform stay close to that of the 
EOB signal up to a frequency fM ^ 0.015, after which 
they start to deviate quite strongly. For reference, note 
that the GW frequency at the 'Schwarzschild-like' LSO 
is equal to M/lso = (67rv^)-i « 0.02166, while the GW 
frequency at the adiabatic-EOB LSO is, for v = 1/4, ap- 
proximately 42% larger, being equal to Mf^°^ w 0.0308. 
In other words, the PN(f) model starts to strongly devi- 
ate from the EOB(f) model at a frequency smaller by a 
factor two than the adiabatic-EOB LSO frequency. Later 
we shall study in more detail after which frequency the 
difference between PN(f) and EOB(f) starts to introduce 
an unacceptable level of 'inaccuracy' (in the technical 
sense defined above). 
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FIG. 2: Upper panel: effective GW amplitudes hs{f) = 
\fh{f)\ (see text) of the various frequency-domain waveform 
models considered in this paper for an equal-mass system 
[assuming a fiducial distance of 1 Gpc and a total mass of 
M = 75Mq, in order to get actual amplitude values]. Lower 
panel: differences between the phases of the closed-form 
frequency-domain waveforms, and the phase of the frequency- 
domain EOB waveform. One has subtracted a linear term 
af + b corresponding to an arbitrary initial time, and an ar- 
bitrary initial phase. The vertical dashed line indicates the 
GW frequency of the "Schwarzschild-like" LSO. The horizon- 
tal segments near the bottom of the two panels indicate the 
shortest frequency intervals over which 70% of the SNR for 
an inspiral signal is accumulated when observed with Adv. 
LIGO (see text). 



Finally, note that the horizontal segments near the 
bottom of the two panels of Fig. [5] indicate (for 
each quoted total mass M) the shortest frequency 
intervals [/iA/,/2M] upon which 70 % of the SNR 
for an inspiral signal {hs(f) oc /^^^^) is accumu- 
lated. Here we used the Adv. LIGO PSD, and 
found that /i = 41 Hz and /2 = 128 Hz span the 
shortest interval such that jj^ {df / f){hs{f)/hn{f))'^ = 
0.7^0°° with Kif) « f-'^'. 



B. Non-spinning phenomenological waveform 
model (PhenVl) 

In this paper we shall consider two different phe- 
nomenological models. [See, however, the concluding 
section for additional results concerning a recently pub- 
lished, third phenomenological model.] The first model 
was built by Ajith et al. (43| by starting from ~' 30 time- 
domain, hybrid non-spinning waveforms in the mass- 



ratio range, 0.16 < v < 0.25. These time-domain hy- 
brid waveforms were constructed by joining together, 
in the time domain, some early 3.5 PN TaylorTl re- 
stricted waveforms with relatively short (~ 4 inspiral cy- 
cles, plus merger) subsequent numerical waveforms gen- 
erated by the Jena group using their bam code [67l.l68j 
with fourth-order finite differencing. Afterwards, the re- 
sulting phenomenological model was tested against some 
other hybrid waveforms built from longer and more ac- 
curate numerical simulations produced using the CCATIE 
[2H and the same BAM code employing sixth-order fi- 
nite differencing. These time-domain hybrid waveforms 
were then Fourier-transformed, and the resulting sam- 
ple of frequency-domain waveforms /ihybrid(/) was then 
fitted to some piecewisc-continuous multi-parameter ana- 
lytic expressions, which define the closed- form 'version- 1 
phenomenological waveform', say 



/iPhenVl (/) = ^PhcnVl {f)e-'^-^^"^^ . (43) 



The phenomenological waveforms arc defined directly 
in the Fourier domain: the amplitude, ^phcnVi(/)j 
is defined as a piecewise function distinguishing inspi- 
ral, merger and ring-down stages, while the phase^, 
i-'PhcnViif), is defined as a PN-type series with phe- 
nomenological coefficients. The whole waveform is char- 
acterized by 10 phenomenological parameters (4 for the 
amplitude and 6 for the phase), each of them fitted to 
a quadratic polynomial of the symmetric mass ratio, v. 
We refer the reader to Sec. IV. C of Ref. [i^ for all the in- 
formation about the actual definition of this phenomeno- 
logical model. 

Since the hybrid waveforms were built using a re- 
stricted PN model, i.e. with a time-domain Newtonian 
amplitude, we expect (in view of the fact that the time- 
domain PN correction is smaller than one during the 
early inspiral) that the amplitudes of the corresponding 
Fourier-domain waveforms will be larger than the ones 
fully incorporating PN effects; and indeed one can see on 
Fig.[2]that ^phcnvi(/) is roughly 6% larger than the EOB 
waveform amplitude. Concerning the phase, note that, 
though full PN information (up to 3.5PN order) was used 
to construct the time-domain hybrid signals, after fitting 
their Fourier transform to closed- form expressions, the 
final frequency-domain phase '(/'PhenVi(/) exhibits signif- 
icant differences with respect to the EOB phase. These 
global discrepancies cannot be absorbed in the arbitrary 
term af -\- b present in the phase, and will be quantified 
below in terms of effectualncss and inaccuracy. 



^ Notice the explicit change of sign of the phenomenological phase 
defined in [4311 compared with our definition here, Eq. I I12I I. 
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C. Spinning, non-precessing phenomenological 
waveform model (PhenV2) 

Recently, Ajith et al. have extended the first phe- 
nomenological model to black hole binary systems 
with non-precessing spins (i.e. spins parallel/anti-parallel 
to the orbital angular momentum) |44{ . Moreover, some 
other improvements have been considered; for instance, 
they are using longer numerical waveforms, covering at 
least eight cycles before merger and also, they arc match- 
ing (in the time domain) NR simulations to more accu- 
rate PN TaylorTl waveforms that include 3.5PN phase 
accuracy, 3PN amplitude corrections to the dominant 
quadrupole mode, and 2.5PN spin-dependent correc- 
tions. The improvement in the amplitude modeling due 
to this new set of phenomenological waveforms is visible 
on Fig. H 

For the purposes of our study, we shall only consider 
non-spinning waveforms, and therefore we shall be us- 
ing just a particular case (x = 0) of the phenomeno- 
logical waveforms defined in [i^l- The definition of the 
closed-form, frequency-domain phenomenological ampli- 
tude and phase, {AphcnV2(/), '0PhonV2(/)}, is very simi- 
lar to what has been explained in the previous section, 
and is given in Eqs. (l)-(3) and Table I of Ref. [44|. 

In what follows, we will refer to this 'version-2' phe- 
nomenological model as PhenV2, in contrast with the for- 
mer PhenVl. The effective amplitude and phase PhenV2 
are compared to FOB, and the other models, in Fig. [5] 



IV. 



FIDUCIAL TARGET WAVEFORM USED IN 
THIS STUDY: NR-CALIBRATED EOB 
WAVEFORM 

A. Effective-One-Body (EOB) formalism 



The FOB formalism [28l - l31j has the unique feature of 
being able to incorporate information coming both from 
PN theory and numerical relativity (NR) simulations into 
an accurate description of the full dynamics, and GW ra- 
diation, of a compact binary coalescence, from the early 
inspiral up to the end of the merger. In addition to in- 
corporating the full PN information available, the FOB 
formalism goes beyond PN theory, even during the in- 
spiral, in several ways: (i) it replaces the PN-expanded 
results by resummed expressions that have been shown to 
be more accurate in various cases; (ii) it includes more 
physical effects, notably those associated with the non 
adiabatic aspects of the inspiral. We recall that all the 
various 'PN models' presently considered in the litera- 
ture make the approximation that the GW phase can be 
computed in an adiabatic manner, by using simply the 
balance between the GW energy flux and the adiabatic 
loss of binding energy of the binary system. On the other 
hand, it has been found in [2^ that the number of GW 
cycles before the LSO over which non-adiabatic effects 
are important is of order ^ {4:iy)~^^^. [Sec Fquations 



(4.47)-(4.49) and Figure 7 in [2i|; and also Ref. ^ for a 
similar result in the v ^ limit.] This indicates that the 
PN phasing will accumulate, because of non-adiabatic 
effects, a dephasing with respect to the FOB one which 
becomes ~ (4z£)~^/^ radians when reaching the LSO (see 
Figure 11 of [231 the = 1/4 case). The accumulated 
dephasing due to the non-resummcd/resummcd differ- 
ence between PN and FOB can easily be larger than this 
non-adiabatic effect. 

Over the last years, it has become possible to improve 
the FOB waveform by comparing it to the results of ac- 
curate NR simulations and by using the natural 'flexi- 
bility' of the FOB formalism to 'calibrate' some FOB 
parameters representing either higher-order perturbative 
effects that have not yet been analytically calculated, or 
non-pcrturbative effects that can only be accessed to by 
NR simulations (33l - [38| . In addition, several theoreti- 
cal improvements have been brought in the FOB formal- 
ism, notably concerning new ways of resumming the GW 
waveform [s^, [zO] , and the G W radiation reaction [s^l • 

The currently most accurate version of the FOB wave- 
form (37j agrees, within the numerical error bars, with the 
currently most accurate NR waveform [l^, [22|, [s^ . This 
means dephasing levels of order of ±0.02 during the en- 
tire inspiral and plunge. In view of the expressions given 
for the 'inaccuracy' functional above, this roughly corre- 
sponds to an inaccuracy level X ~ (0.02)^ = 4 x 10~*. As 
we shall find below that the inaccuracy, with respect to 
FOB, of the closed- form models presented in the previ- 
ous section is much larger than this level, we can mean- 
ingfully use the continuous bank of FOB waveforms as 
'fiducially exact' target waveforms in order to gauge the 
cffcctualncss, and the accuracy of the closed-form mod- 
els. More precisely, we shall use the FOB model de- 
fined in Ref. [13] with FOB parameters 05 = —6.37 and 
flg ~ 50 and with a Pade (3,2) resummation of the func- 
tion f22{x). 



B. Time-domain EOB waveform 

From the practical point of view, the FOB waveform 
is defined in the time domain by integrating some ODF's 
that cover the inspiral down to the 'FOB merger', and 
then by attaching a ring-down signal, in a smooth man- 
ner, at the instant of the 'FOB merger'. [Note that the 
time variable used in the FOB model is the dimension- 
less variable t = t/M. The only parameter entering the 
integrated ODF's is the symmetric mass ratio v.] The 
FOB model is therefore 'semi-analytical' in the sense that 
it needs to use some numerical integration of the basic 
analytical ODF's defining the FOB dynamics. We have 
used a Matlab implementation of the FOB dynamics that 
has not been optimized for speed and/or efficiency. The 
integration time is typically less than 0.15 s per orbital 
cycle [value obtained running on a single core of a 2.80 
GHz quad-core Intel Xeon desktop machine with 4 Gb of 
RAM] , so that we are able to simulate several thousands 



13 



of orbital cycles in less than half an hour by running on 
a single CPU of a desktop machine^. Given its accu- 
racy and speed, the EOB model is a very good candidate 
to be used as a full (inspiral-merger-ringdown) waveform 
generator in actual searches of GW, especially as it is 
clear that some optimization would allow to significantly 
reduce the current integration time per orbital cycle. 

Though this procedure defines extremely long time- 
domain EOB waveforms, it delivers them with an abrupt 
beginning, say at the (dimensionless EOB) time 4 (= 0). 
This abrupt beginning at t}, is a necessary consequence 
of having started the EOB evolution at some finite ini- 
tial radius tq. [Note in passing that, though we start 
the EOB waveform with some specified initial phase at 
^EOB = 0, the fact that we then continuously vary param- 
eters such as ro and v effectively implies that the phase 
of the generated EOB waveform around coalescence is 
essentially random.] 

A direct Fourier transform of the (time-windowed) 
EOB waveform would introduce, because of edge effects, 
some significant oscillations in the frequency domain that 
would "spoil" the beginning of the frequency-domain 
waveform. [See, e.g. , 0| for a detailed discussion of 
the oscihations in h{f) induced by a sharp edge in h{t).] 
In order to lose as little information as possible from the 
original EOB evolution, we have "smoothed" the initial 
edge of the time-domain waveform coming out of the 
numerical integration by extending it to earlier times, 
t < t},. Namely: 

1. Since the times t < it arc deep into the adiabatic 
inspiral, the signal there would be analytically well 
described by some post-Newtonian chirping time- 
domain model. However, as our aim is purely tech- 
nical, namely to avoid the oscillations due to edge 
effects to spill over frequencies / larger than the 
initial time-domain frequency F{tb), we found suf- 
ficient to move the abrupt starting edge towards 
earlier times by attaching, say for ta < t < tb, a 
simple phenomenological Newtonian time-domain 
evolution of the form 

A{t) ^ a2{ai^i)-^/^ 

We determined the 5 phenomenological parameters 
an,§n by imposing continuity, aX t = of A, A, 
(j), 4> and (j) with the EOB. Given the one-to-one 
relation between frequency and time in a chirping 
signal, this oversimplified left extension is enough 
to displace the edge effects due to the abrupt start- 
ing of the waveform to a controllable frequency- 
domain region localized around F(ta). We choose 



By comparison, note that a typical NR simulation takes about 
two days per orbital cycle on a computer cluster with several 
hundreds of cores, which is a factor ~ 10* longer [in actual CPU 
time] than EOB. 



ta SO that F{ta) is sufficiently smaller than F{tb), 
thereby having a smooth frequency-domain wave- 
form for any / > F{tb). 

2. Finally, to further suppress any edge effect related 
to the sharp start of the waveform at the left- ward 
displaced time ta, we multiply this Newtonian-like 
time-domain waveform by a (time-domain) semi- 

Tukey filter in the interval ia < i < ia + 3*° . 
This filter replaces the sharp start at t ~ ta by 
a progressive build up of the waveform during the 
first third of the Newtonian extension. 

At the end of this procedure, we have a very long, 
smooth time-domain waveform that vanishes for t < ta, 
smoothly starts building up at i = ia, and ends up de- 
caying exponentially [deep in the ring-down stage] when 
t +00. This time-domain waveform is now ready to be 
Fourier transformed. [A similar procedure has been used 
by authors who needed to compute the Fourier transform 
of short, time- windowed time-domain NR waveforms, see 
e.g. m.] 



C. Frequency-domain EOB waveform 

Starting from the extended time-domain EOB wave- 
form defined in the previous subsection, we now compute 
its EFT. However, we must discard from the resulting 
frequency-domain waveform the low-frequency contribu- 
tion coming from the Newtonian extension, i.e. anything 
from fM < fbM, where fbM = 0EOB(4)/(27r). This 
means that, so far, we have constructed only a version of 
the full frequency-domain EOB waveform which is cut-off 
on the low-frequency side: /ieob(/)^(/— /&)■ We wish now 
to extend to lower frequencies, f < fb, this frequency- 
truncated EOB waveform in order to have the complete, 
frequency-domain EOB waveform. We can do this low- 
frequency extension by joining together, in the frequency 
domain, an early frequency-domain FN waveform PN(f) 
(as defined above) to the frequency-truncated EOB wave- 
form /ieob(/)^(/ — fb) we just constructed. We will do 
this joining of PN(f) with EOB(f) in a smooth and pro- 
gressive manner, over a certain finite frequency interval 
[/o,/i], entirely located within the range of frequencies 
where we have, so far, computed /ieob(/)- In other words, 
we impose the inequalities /b < /o < /i- The choice 
fa ~ fb will define the best approximation we can con- 
struct of the complete frequency-domain EOB waveform 
EOB(f). 

Technically, we accomplished this progressive 'joining' 
between a low-frequency PN(f) and a higher-frequency 
EOB(f) by constructing 

HfM) = (1 - e/„M,/,M(/M)) hl^^^^ifM) 

+ QhM.,hM{IM)h^o^{fM) (44) 
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where Q faMjiAiix) is an infinitely smooth, C°° , step- 
like transition function that goes from to 1 in a finite® 
interval [/qM, /iM]. Specifically, we have used 



r 



1 



X < Xq 

xq < X < xi (45) 

X > Xi 



In addition, we have chosen the extrinsic parameters 
tc, 0c entering the PN(f) waveform hpn{f) (see Eq. (|40|) ) 
so as to minimize the phase difference between ipPNif) 
and '!/'eob(/) in the junction interval. Note that the PN(f) 
waveform, hpf,{f) is generated with the same physical pa- 
rameters {DcS, M, i^} that were used for the EOB simu- 
lation. 

Finally, for the practical purpose of this study, we shall 
define the complete frequency-domain EOB waveform 
EOB(f) as the 'smooth frequency-domain PNU/p j^EOB 
hybrid' constructed in Eq. pi)) , with the choices /o ~ /t 
and /i = 1.05/o, where /{, is the GW frequency corre- 
sponding to the beginning of the longest EOB simula- 
tions we generated. More precisely, we have explored 
a fine grid of mass ratios q = 7712 /mi > 1 between 
q ~ 1 and q ~ 10. Given our current implementation 
and running our simulations on a single desktop com- 
puter with a 2.80 GHz quad-core and 4 Gb of RAM, 
we could, when I < q < i, start the simulations at 
To = Ro/M = 90, which corresponds to an initial GW 
frequency Mfb = l/^irrQ^^) ss 3.73 x 10~^. For larger 
mass ratios, up to (7 = 10, we started at vq = 76, which 
corresponds to Mfb « 4.80 x 10^*. The effective ampli- 
tude and phase of this complete EOB(f) waveform is rep- 
resented in Fig.[3]for the values q = 1,2, 4, 10 of the mass 
ratio. Note that this represents simulating 1553, 1744, 
2414 and 3042 orbital cycles, respectively, that took be- 
tween 4 and 10 minutes of computational time on a single 
CPU. 

The starting point of these simulations is in the very 
deep inspiral phase of the evolution, where both SPA and 
PN approximations are perfectly valid. Thus, we shall 
consider these waveforms as the "correct" ones, using 
them as the target when computing the inaccuracy of 
PN-SPA. 



V. ACCURACY UPPER LIMIT ON THE 
JUNCTION FREQUENCY /„ OF HYBRID 

PN(/) U/„ EOB(f) WAVEFORMS 

In this section we shall address a question that is a 
side-issue with respect to the main aim of this paper. 
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FIG. 3: EOB(f): frequency-domain gravitational wave eflec- 
tive amplitudes, ,f\h{f)\ (top panel), and phases, V'(./) (bot- 
tom panel) , obtained from EOB simulations of BBHs coa- 
lescences with different mass ratios. [As in Fig. [21 the actual 
values of the amplitude have been obtained assuming a total 
mass of 75Mq at 1 Gpc] 



but which is important in the general context of con- 
structing accurate, frequency-domain waveform models. 
This question concerns the upper limit on the frequency 
/o parametrizing^ the joining between a lower-frequency 
PN(f) waveform to a higher-frequency EOB(f) waveform 
which does not lead to an undue loss of accuracy for the 
resulting hybrid waveform PN(f)U/QEOB(f). We shaU 
measure the loss of accuracy of the (frequency-domain) 
hybrid PN(f)U/jjEOB(f) by means of the 'inaccuracy 
functional' T[hm[fo]', h^], Eq. (P^)) . where /im[/o] denotes 
the hybrid PN(f)U/oEOB(f); and where denotes our 
best approximation to the full, exact EOB(f) waveform, 
i.e., as explained above, the hybrid PN(f)U/^EOB(f) gen- 
erated with the very low frequency fb at which we started 
the EOB simulation. We recall that Mfb ~ 4 x IQ-" cor- 
responds to such a non-relativistic and adiabatic part 
of the inspiral (thousands of orbital cycles before the 
merger of the BBH) that we can safely consider the lower- 
frequency PN(f) piece in the hybrid PN(f)U/,EOB(f) as 
being indistinguishable from the low-frequency part of 
the exact EOB(f). On the other hand, in view of Fig. [21 
we expect that PN(f) will start to significantly deviate 
from EOB(f) for frequencies somewhere between fM ~ 
10"^ and fM lO'^, thereby leading to an increasing 
inaccuracy I[h^[fo]; h^] of the hybrid PN(f)U/oEOB(f) 



Note that the often used smooth transition functions based on 
hyperboHc tangent functions are also C°° , but interpolate be- 
tween to 1 only over the infinite real line. 



^ For simplicity and concreteness, we shall henceforth use, for the 
smooth joining of PN(f) with EOB(f), Eq. 1 144 I I. a frequency 
interval [/o,/i] with a fixed ratio /1//0 = 1.05. This leaves 
the single parameter fo to characterize the 'junction frequency' 
between the two waveforms. 
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as foM reaches such values. 

There are two motivations for this study. First, though 
we have been able, for this work, to easily generate many, 
very long time-domain EOB waveforms, containing thou- 
sands of orbital cycles, and then use an automated ver- 
sion of the procedure explained above to generate the cor- 
responding frequency-domain EOB waveforms, it might 
be useful, when one will need to construct dense banks 
of accurate frequency-domain EOB waveforms, to cut 
down to a strict minimum the integration time of each 
EOB simulation. Determining the maximum /o lead- 
ing to some given inaccuracy will be useful in allowing 
one to start the EOB simulation at a higher GW fre- 
quency /{, = /o- A second motivation, is that our result 
will clearly also tell us the maximum frequency at which 
we must start a NR simulation {i.e. the minimum num- 
ber of inspiral orbital cycles it must include) if we wish 
to be able to extend (to lower frequencies) its Fourier- 
transformed GW signal by using a simple closed-form 
PN(f) waveform. 

As said above, in order to study the frequency up- 
per limit for the accuracy of PN(f) waveforms, we gen- 
erate frequency-domain hybrid waveforms h,n[fo]{f) ~ 
PN(/) U/j, EOB(/), and compute the relative inaccuracy 
I[hrn[fo]'i hj-n[fb\] ^ fuuctiou of the joining frequency 
Mfo, or, equivalently, AIluq = 27rAf/o. The inaccuracy 
functional will also depend on the particular symmet- 
ric mass ratio v one is considering, and on the value of 
the total mass M (in seconds), since the detector noise 
curve that enters the definition Eq. ([2^ depends on the 
physical frequency / (in Hz) and not on the dimension- 
less scaled frequency fM. In particular, we shall consider 
four different values of h' corresponding to the mass ratios 
q ~ 1,2,4,10 and we shall maximize I{Mfo,i/,M), for 
each value of Mfo, and v, over M, i.e. we shall consider 
the worst case scenario for the loss of accuracy. Since we 
arc considering ground-based detectors, the mass range 
wiU be considered to be M € [3, 5OO]M0. This yields, for 
each value of the worst-case inaccuracy as a function 
of M/q. We have done this study for the three different 
detector noise curves discussed above, and represented 
in Fig. [T] initial LIGO, advanced LIGO, and advanced 
Virgo. Note that the inaccuracy functional only depends 
on the shape of these noise curves, and not on their ab- 
solute magnitudes. 

We represent in Fig. U] the inaccuracy, 
^[hm[fo]]h„i[ff,]], maximized over M, as a function 
of the junction frequency Mfo^ for 4 different mass 
ratios, and for the Advanced LIGO noise PSD. Accord- 
ing to Eq. pop , the maximum inaccuracy that we can 
tolerate depends both on the signal-to-noise ratio, p, and 
on the choice of an additional "safety factor" e. More 
precisely, what matters is only what we can call the 
effective SNR, namely poff = p/e- Two models will be 
indistinguishable if their relative inaccuracy is smaller 
than 1/Poff- present section we shall consider the 

effective SNRs peS = 10 or p^s = 20. Note that, as we 
expect typical observations to have a real SNR of 10, 



the inaccuracy limit l/p^ff with p^g = 10 corresponds to 
taking e = 1, i.e. to allowing the inaccuracy of the model 
waveform to double the effect of the detector noise. It 
would seem more reasonable to take at least e = 1/2, 
i.e. to consider an effective SNR of pcff = 20. 

The most interesting feature of Fig. U is the notable 
difference between the case of mass ratios near unity 
(namely q = 1 and g = 2), and the case of large mass 
ratios: q = A and q = 10. When g = 1,2 Fig. |4] shows 
that the accuracy loss remains quite small up to frequen- 
cies rather close to the Schwarzschild-likc LSO frequency 
M/lso = (67rv^)-i w 0.02166. [One should, however, 
recall that the real (EOB-predicted) LSO frequency is 
42% larger, being close to M f^°^ w 0.0308.] In other 
words, the PN(f) model produces accurate waveforms up 
to near the LSO when q = 1,2. By contrast, when the 
mass ratio increases, namely 9 = 4, 10, the inaccuracy 
gets significant for frequencies which are much smaller 
than the LSO frequency (especially if one considers the 
threshold pcfi = 20). As the number of orbital cycles be- 
tween /o and /i is approximatively given by the "New- 
tonian" estimate (see, e.g. , [soj ') 

A^^'^M/o,/!] = [(7rM/o)-^/3 - (7rM/i)-5/3]/(647rz.) 

(46) 

such a smaller M/o means a much larger number of or- 
bital cycles left between /o and, say, /i ^ /lso (both 
because of the power-law dependence on Mfo, and of 
the explicit proportionality to 1/j^). Even if we consider 
the somewhat lax threshold pcS = 10, we see on Fig. |3] 
that the PN(f) waveform can be used: up to Mf « 0.014 
when g = 1 or 2; only up to a twice smaller frequency 
Mf « 0.008 when (? = 4; and only up to a seven times 
smaller frequency Mf « 0.002 when q = 10. 

The Table of Fig. 2] gives more quantitative informa- 
tion about this accuracy upper limit on the junction fre- 
quency to PN(f). The table considers all three detectors 
noise PSD, and the two thresholds peff = 10 and peS ~ 20 
(the quantitative data pertaining to the pcs = 20 thresh- 
old are given within parentheses). In particular, the last 
two columns indicate the conversion of M fo in separa- 
tion tq ~ Ro/M between the BHs. and in number of 
orbital cycles left before reaching the merger. [The con- 
version between Mfo and rg = Ro/M was formally done 
using the approximate formula r = (tt/M)"^/"^, while 
the number of orbital cycles was computed by using the 
underlying EOB dynamics; defining the "merger" as the 
instant where the (metric) waveform reaches its maxi- 
mum amplitude.] 

Notice the steep rise of the number of orbital cycles left 
as the mass ratio gets larger than 4. For a given accuracy 
threshold 1/Peff the Advanced Virgo detector is the most 
demanding one for the PN(f) model because the shape 
of its effective noise curve peaks at a lower frequency, 
but note that, for a given signal amplitude the SNR of 
Advanced Virgo would be smaller than that of Advanced 
LIGO. 

Leaving to future work [7l[ a detailed discussion of the 
precise meaning of these results (notably in terms of the 
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intrinsic accuracy of the PN(f) waveform), let us men- 
tion that our resuhs suggest that the "inaccuracy" of a 
frequency-domain hybrid, say PN(f)U/oNR(f), between 
a lower-frequency PN(f) and a higher-frequency NR(f) 
waveform will be similar to the one plotted in Fig. |4l In 
particular, we think that, if one requires the same type 
of threshold on the inaccuracy, the number of orbital cy- 
cles before merger will have to be similar to the numbers 
indicated in the last columns of the Table. On the other 
hand, let us recall that numerical simulations arc compu- 
tationally very expensive, and the current longest wave- 
forms being produced span at most the last 10—15 orbital 
cycles of the evolution (tI, [t^. Recently, a joint effort 
between several NR groups and scientists producing an- 
alytical models [t^I has begun to make use of TeraGrid 
resources in order to produce longer and more accurate 
simulations. However, one docs not currently foresee the 
possibility of producing more than 15 — 20 orbital cycles. 
Comparing these current NR capabilities with the num- 
bers we are obtaining here for the upper-frequency-limit 
of frequency-domain hybrids using PN(f) for / < /o, we 
observe that, for mass ratios g > 4 it will be impossi- 
ble to "join" a PN(f) waveform to a NR(f) one with- 
out incurring an unacceptably large accuracy loss. In 
intuitive terms, one can say that we are here facing the 
probable existence of a frequency gap where neither PN 
theory nor current NR simulations are able to provide ac- 
curate enough information. We think, however, that the 
effective-one-body approach, aided by suitably designed 
NR calibrations, may be able to bridge this dangerous fre- 
quency gap, and then to construct some accurate-enough 
complete waveform models. 

VI. VALIDITY RANGES OF PN(/) AND OF 
PHENOMENOLOGICAL MODELS 

This section will present the central results of our 
work: namely the quantitative study of the ineffectual- 
ness and of the inaccuracy of the closed-form, frequency- 
domain waveforms presented above, with respect to the 
fiducial target waveform defined by the NR-calibrated 
EOB waveform (computed in the frequency-domain, as 
explained above). [For brevity, we shall sometimes de- 
note the NR-calibrated EOB waveform as EOBnr-] Our 
aim is to assess the regions of the parameter space of 
(non- spinning) BBHs, i.e. M and v, where: 

• either one of the closed- form models (PN(f), 
PhenVl or PhenV2) provides an effectual represen- 
tation of EOBnr; as measured by the value of the 
effectualness, see Eq. (jl9p : or, 

• either one of the closed- form models (PN(f), 
PhenVl or PhenV2) provides an accurate represen- 
tation of EOBnr; as measured by the value of the 
inaccuracy., see Eq. (|30[) . 

We did not investigate whether a suitable 
PN(f)U /oPhenV2(f) hybrid might be closer to EOBNR(f). 



Indeed, Phen waveforms are already supposed to have 
incorporated a PN-type early behavior, with an NR-type 
later waveform. 

Before presenting our quantitative results, let us an- 
ticipate some general aspects of our results by noting 
that one a priori expects that neither PN(f) nor phe- 
nomenological models will be able to accurately describe 
all stages of the evolution of BBHs. Through the detec- 
tor noise PSD, this lack of uniform physical validity will 
then translate into validity ranges in the space of physical 
parameters, M and ly. In particular, we expect the PN(f) 
waveforms to only provide a valid description of the early 
inspiral stage of the evolution up to a certain frequency 
[~ fo{M, v), see Sec. IVl above], but to start failing as the 
two BHs get close and the dynamics becomes more rel- 
ativistic and less adiabatic. Given the inversely propor- 
tional relation between the total mass and the frequency, 
inspiral-only models, such as PN(f) (or, for that mat- 
ter, any other PN waveform), will be adequate only for 
searches of relatively low mass systems, M S [2,2O]M0, 
in initial ground-based detectors (for such systems the 
merger occurs at higher frequencies than the detectors 
sensitivity window; see the horizontal segments near the 
bottom of Fig. [5]) . On the other hand, the phenomeno- 
logical models are built to effectively describe the merger 
phase of the evolution. They can therefore be expected 
to provide an adequate method to search for high mass 
systems (say M > 50Mq), where the merger phase is the 
main contribution to the SNR (see Fig. [5]). This prelimi- 
nary discussion suggests that there might exist an inter- 
mediate range of masses, around 20Mq ^ M < 5OM0, 
where none of the closed-form models can provide an ac- 
curate description. We shall sec that this is indeed the 
case. 

Our quantitative results on the validity range of these 
two models, both from the effectualness (detection) and 
accuracy (measurement) points of view (using EOBnr 
waveforms as fiducial targets), are presented in Figs. [5] 
and [HI and in Table HI We now explain them in detail. 



A. Effectualness. Detection 

As recalled above, when one is only interested in de- 
tecting GW signals from compact coalescing objects but 
not in the extraction of accurate physical parameters 
from the observations, an important indicator to eval- 
uate a given approximated model is the effectualness, £, 
Eq. ((Ti)) , which represents the fractional loss of SNR en- 
tailed by the use of a sub-optimal match filter. Actually, 
the most important indicator is the "effective effectual- 
ness" which takes into account the discrete character of 
the template bank (see above). We have indicated in 
Eq. the range of thresholds that one usually consid- 
ers for the effectualness. Let us say again that the often 
quoted effectualness level of 97% (corresponding to an 
incffectualness £ = ! — £< 3%) is a very minimal re- 
quirement that does not seriously take into account the 
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FIG. 4: (Plot) "Inaccuracy" l[hm[fo]; hm[fb]] [defined by Eq. ^ or Eq. ^] of frequency-domain hybrid PN(f)U/oEOB(f) 
waveforms as we move their junction frequency, Mfo (= Ma;o/27r), to higher values; taking as reference waveform an hybrid 
PN(f)U/(,EOB(f) where the junction frequency is extremely low {Mfb ~ 4 x 10"*). We have computed I for 4 diflerent 
mass ratios, using the Adv. LIGO noise PSD and maximizing the inaccuracy over M £ [3, 5OO]M0. The vertical dashed 
line represents the frequency of the Schwarzschild LSO and the two horizontal dashed lines indicate the accuracy thresholds 
corresponding to the indicated effective SNRs pcff = p/e, according to Eq. H30p . (Table) Upper limits (for three different 
detectors noise PSDs, and four mass ratios) on the value of the junction frequency Mfo above which the inaccuracy of the 
hybrid PN(f)U/QEOB(f) becomes larger than the threshold corresponding to an effective SNR equal to 10 and 20. We have 
also converted them into a minimum distance between the BHs, Ro, and a minimum number of orbital cycles before merger 
(see text). 



constraints coming from the discrete nature of the tem- 
plate bank. It would seem more appropriate to require 
an ineffcctualness near £ < 1%, which, combined with an 
equivalent minimal match smm = 1%, guarantees that 
one does not lose more than 6% of the potential events. 

To have a first view of the influence of choosing an 
effectualness (and accuracy) threshold on the validity 
range of various models, the top panels of Fig. [5] plot 
as dark gray (blue online) lines the ineffectualncss £ as a 
function of the total mass, for a fixed mass ratio of 4:1, 
and using the Adv. LIGO noise PSD. The three top pan- 
els respectively correspond to the three closed-form mod- 
els h„i introduced above: PN(f), PhenVl and PhenV2. 
[We recall that EOBnr is used as fiducial "exact" target 
hx-] 

These top panels show that as one imposes a stricter 
threshold on the effectualness {i.e. a smaller level for the 
ineffectualness £ = ^ — £) the validity range of the var- 
ious models change in the direction expected from the 
preliminary discussion above. For the PN(f) model (left 
top panel), the validity range at some ineffectualness 
threshold £ is an interval of the type M < M™^^{£), 
where the maximum allowed total mass M^^^{£) de- 
creases as one decreases £. For instance, at the level 
£ = 0.03, indicated by the upper dashed line in the left 



top panel (blue online), one has M™^^ 17Af0, while 
the stricter level £ = 0.01 leads to M™^^ sa UMq. This 
behavior is consistent with the fact that the PN model 
is only adequate during the inspiral, up to some max- 
imum dimensionless frequency / = fM < fQ{v) (see 
Fig. U). Then, given the fixed detector sensitivity win- 
dow in Hz, this translates into an upper limit on the 
total mass: M < /o('^)//dotcctor- On the other hand, the 
phenomenological models considered in this paper were 
built with the specific aim of reproducing the late inspi- 
ral, merger and ringdown phases of the evolution found 
by NR simulations, and therefore, they should be most 
effectual for the detection of high mass systems. And, 
indeed, one finds that, for the phenomenological models 
(right top panels) the validity range at some ineffectual- 
ness threshold £ is an interval of the type M > M^°''{£), 
where the minimum allowed total mass M^^"^{£) in- 
creases as one decreases £. For instance, at the level 
£ = 0.03, and for PhenV2, as indicated by the upper 
dashed line in the rightmost top panel (blue online), one 
has ^ UMq, while the stricter level £ = 0.01 

leads to a much larger value: M^f^""^^ sa 41Mq. [For 
PhenVl, these minimum masses are M^f"^^^ « 16Mq 
for £ = 0.03; and M^^^^i ^ 246Mq for £ = 0.01]. This 
example (with q = A and Adv. LIGO) shows that: (i) 
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FIG. 5: Top panels: InefFectualness £ — 1 — £ [dark gray (blue) solid line; see Eq. (|14p ] and inaccuracy I [light gray (red) 
solid line; see Eq. (|27|) or Eq. (|29|) ] of the different closed-form models compared to an EOB(f) waveform, as functions of the 
total mass M of the system (for a 4:1 mass ratio, and using Adv. LIGO noise PSD). The horizontal dashed lines represent 
rather mmimal thresholds on these quantities: 3% for the ineffectualness, and 1% for the inaccuracy (see text). Bottom panels: 
Fractional biases on the M and v parameters after maximizing the overlap between the model and target waveforms when 
computing the effectualness (see Eq. (|15[) '). The horizontal dashed line represents the zero biases and the vertical one denotes 
the M value where the ineffectualness reaches 3% in the top panels. 



at the minimal effectualness level E = 0.03 there is an 
overlap between the 'low-mass' interval where PN(f) is 
effectual, and the 'high-mass' interval where PhenV2 is 
effectual, so that the union of the two template banks 
is able to detect the full mass range; while, by con- 
strast, (ii) at the stricter effectualness level £ ~ 0.01, 
there exists a mass gap, between M^^^ « 14M0 and 
^Phe„V2 ^ 41M0, in which neither PN(f) nor PhenV2 
(nor Phen VI) is effectual enough. We shall focus on 
these mass gaps below, after having commented upon 
the other panels of Fig. [SJ 

The bottom panels of Fig. [S] complete the discussion 
of the effectualness calculation by showing the fractional 
biases between the actual M and v parameters and the 
ones that maximize the overlap between the model and 
the target waveforms in the maximization, Eq. (jl5|) . In 
our analysis, we have restricted the parameters to have 
a physical meaning, i.e. v < 0.25. However, the results 
obtained by allowing the parameters to be unphysical re- 
main almost the same: only a small improvement of the 
effectualness near the equal-mass region is observed, be- 
cause this is where the restriction becomes relevant. By 
looking at our particular results in Fig. [5] for v = 0.16 
and the Adv. LIGO noise PSD, we observe that the sys- 



tematic biases in estimating the parameters are smaller 
than 1% for M < 25Mq when using PN waveforms; 
while they are smaller than 5% for M (10% for v) for 
AI e [8, 2OO]M0 when using phenomcnological models. 
Anyway, these fractional biases are merely informative 
since one would have to compare them with typical sta- 
tistical errors from a proper parameter estimation study 
in order to truly interpret the consequences of having an 
effectual, but not accurate, model. Let us mention that 
several works in the literature have studied the param- 
eter estimation of BBH coalescences with ground-based 
detectors, using cither PN templates [s^, [75l - [80| or 
phenomcnological ones [56j . 

Let us now study in detail the mass gaps where 
none of the closed-form models arc effectual enough. 
Since the non-preccssing phenomcnological waveforms 
(PhenV2) [i^l are always more accurate and produce 
smaller mismatches with FOB than the original non- 
spinning ones (PhenVl) (43j . we shall study, for various 
effectualness levels, the region of parameter space where 
the union of PN(f) and PhenV2 models fails to be effec- 
tual. As we saw above, this region is defined, for each 
given ineffectualness level £, and for each given value of 
ly, by the mass gap between M™^^ and Mj^jj™^^, i.e. by 
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TABLE I: Total mass ranges where the ineffectualness (Eq. (|14[) ') and inaccuracy (Eq. (|30|l 'l of both PN(f) and PhenV2 models 
are greater than the quoted numbers in the top row of each table. The interpretation of these ranges are the following: (Left) 
Mass ranges where both models have an ineffectualness larger than the quoted numbers; (Right) Mass ranges where both 
models have an inaccuracy larger than 1/Pcffi where the effective SNR paS — {5,7,10,20} respectively. An asterisk in front 
of a mass range denotes that a small range in between the quoted M values does not satisfy the quoted inequalities for the 
ineffectualness / inaccuracy. 



the interval Af™^(f,i/) < M < M^f^"^^ {£ , ly) , when it 
is non empty. These mass gaps are listed in Table 1(a) as 
three sub-tables for the different detectors: initial LIGO, 
advanced LIGO and advanced Virgo; and in there the 
intervals [M^^^{£ , ly) , M^^^""^^ {£ , v)] are given for four 
different mass ratios (q = 1,2,4, 10) and for four differ- 
ent effectualness levels: 3%, 2%, 1%, 0.45%. 

In Table 1(a), one can observe how these "ineffec- 
tualness intervals" [M^^^^^{£ , , M^^l'^"^^ {£ , ly)] become 
wider as we go to more restrictive effectualness bounds, 
and also as we consider more sensitive detectors. Con- 
cerning their dependence upon v, we see that mass ratios 
between 2:1 and 4:1 generally lead to the narrowest in- 
tervals, while the other mass ratios (especially q = 10) 
lead to wider intervals, mainly because of the increase 
of their right end: M^f^'^''^^{£ , v). Typically, the lack of 
effectualness of the closed-form models takes place for in- 
termediate mass systems, from a total mass of IQMq to 
20 — IOOM0, depending of the mass ratio. For instance, 
considering the Adv. Virgo noise PSD and allowing a 
maximum loss of SNR of 2%, the analytical models would 
not be valid to search for equal-mass systems between 
(4.5 -h 4.5)M0 and (26 + 26)Mq; for 2:1 systems from 
(3.5 + 7)Mq to (5 + 10)Mq; for 4:1 from (2.4 4- 9.6)Afo 
to (4 -I- 16)Mq and for the most extreme mass ratio 10:1 
from (1.6 + 16.4)Mo to (14 + 141)Mq. 

Finally, in order to have a better view of the regions 



of the {M, v) parameter space where the union of PN(f ) 
and PhenV2 fails to be "5-effectual" {i.e. where the in- 
fimum of the ineffectualnesses of PN(f) and PhenV2 is 
larger than £), we have determined the ineffectualness 
intervals [M™^(f , v) , M^f^^"^^ {£ , i/)] on a finer grid res- 
olution in the whole parameter space. The contour plots 
in the large top panels of Fig. [S] and Fig. [7] represent the 
result of this ineffectualness calculation when one con- 
siders the two closed-form models simultaneously. By 
contrast, the small sub-panels located below the large 
upper panels represent the separate values of the effec- 
tualnesses of the two closed-form models: PN(f) on the 
left and PhenV2 on the right. The four contour lines 
represent the same effectualness thresholds used in Ta- 
ble U (i.e. 3%, 2%, 1%, 0.45%), but the higher resolution 
in v allows us now to observe the continuous evolution 
of the validity ranges as the mass ratio varies. As be- 
fore, normally the left-most bounds are determined by 
the total mass values M™^^{£, v) above which the PN(f) 
model is not valid anymore, whereas the right-most lim- 
its are given by the minimum mass values M^™^^(£ , v) 
where phenomenological models are effectual. Choosing 
a certain ineffectualness level £, the corresponding re- 
gions (from the darkest one in which £ > 3% to the 
lightest [but not white] one where £ > 0.45%) on the 
contour plots represent the domain of parameter space 
where none of these closed-form models could be used 
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FIG. 6: Contour plots of the same mass ranges as in Tabled but plotted with a finer mass ratio grid between 1:1 and 10:1 
and also considering the three detector noise PSDs. The top panels plot the regions of the parameter space where the infimum 
of the ineffectualnesses of the PN(f) and PhenV2 models is larger than {3%, 2%, 1%, 0.45%} going from the darkest to the 
lightest areas; whereas the bottom panels show the parts of the parameter space where the infimum of the inaccuracies of the 
PN(f) and PhenV2 models is larger than 1/Poff i where the effective SNR is potr = {5, 7, 10, 20}, when going from the darkest to 
the lightest region. The small sub-panels located below the larger ones represent the results obtained when considering each 
closed-form model separately (with PN(f) on the left and PhenV2 on the right). Here we only consider BBH systems with 
individual components heavier than I.5M0; this excludes the dashed area in the lowest mass range of each plot. The four 
parallel dashed lines represent the locus of the systems where (at least) iVo.gg = {100, 50, 20, 10} orbital cycles (going from left 
to right) contribute 99% of the SNR (see text). 



in an actual search, even for the mere purpose of detec- 
tion^ without requesting a faithful recovery of the physi- 
cal parameters. Let us start by noting that, in the case 
of initial LIGO, the simultaneous use of PN and PhenV2 
waveforms (i.e. , in fact, the use, for each particular Af, v 
of the best of the two, namely the one having the small- 
est inefFectualness for these values of M and ly) allows 
one to cover the full parameter space with an inefFectu- 



alness < 3%, i.e. an effectualness larger than 97%. At 
higher effectualness levels, like 98% or 99%, one is pro- 
gressively losing the complete coverage of the parameter 
space. It is, however, interesting to see, on the smaller 
sub-panels below the top left panel, the separate perfor- 
mances of PN and PhenV2 waveforms: PN has a good 
effectualness only in the small-mass corner on the left. 
On the other hand, PhenV2 (for initial detectors) is 99% 
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FIG. 7: Same plots as top panels in Fig. [6] (ineffectualness), but considering PN(f) - PlienVl models instead. In this case, 
the accuracy results have not been plotted because neither PN(f) nor PhenVl are accurate for M > 20Mq at any mass ratio 
(between 1:1 and 10:1), even when considering an effective SNR of 5. 



effectual in a large part of parameter space, though it 
fails to be 98% effectual in three sub-zones: (1) in the 
small-mass corner, (2) in an "island" centered around 
M — IOMq, q = 1, and (3) in another island located 
above g > 9. This good effectualness performance of the 
closed-form models is, however, drastically decreased in 
the case of advanced detectors. Note, for instance, how, 
in the case of advanced LIGO, the PhenV2 waveform 
fails to be 97% effectual in an extended "C-shaped" re- 
gion which covers many physically interesting systems: 
notably the range 3Mq < M < 40Mo,(j 1, a do- 
main located in the high q region, as well as essentially 
all systems with M < IOA/q. However, the simulta- 
neous use of PN and PhcnV2 waveforms allows one to 
substantially reduce the extension of this ineffectualness 
domain. As seen on the middle top panel, the union of 
the two waveforms fails to be 98% effectual only in two 
"islands" located either around M ~ 20MQ,q ~ 1, or 
around M ~ AOMq, q > 7. Note also, how, in the case of 
the advanced Virgo detector, requiring an ineffectualness 
< 2% {i.e. an effectualness > 98%) excludes the joint use 
of closed-form waveforms in an extended "ridge" which 
separates a low-mass region on its left from a high-mass 
region on its right. 

Before discussing the "inaccuracy analogs" of the inef- 
fectualness regions plotted in the top panels of Fig. [6l let 
us mention the meaning of our results for the most recent 
LIGO- Virgo searches for low-mass compact binary coa- 
lescences (CBC) using PN templates. In particular, I60l- 
[6^ have performed searches for M G [2,35]M0 and |63l| 
for M G [2,4O]M0. Upper limits for the particular sys- 
tems (1.35 + 1.35)Mo, (5 + 5)Mo and (5 + 1.35)M0 
have also been set in [6l|, HI]. Our results show that 



the PN templates used in all the searches done so far 
(with the initial LIGO noise PSD) had an ineffectual- 
ness < 1% for systems with M < IOMq, and that 
this performance could be extended to most of the ana- 
lyzed parameter space by also including phenomcnologi- 
cal waveforms. However, in a physically significant part 
of the parameter space (namely systems with a total mass 
M > IOA/q and mass ratios, either 9 < 2 or g > 5), 
the ineffectualness of the PN templates exceed £ = 1%, 
and even £ = 2%. Keeping in mind that the discretized 
template bank used in the searches had a maximal mis- 
match of Emm = 3%, this means that, in a significant 
part of parameter space, the searches could be missing 
l-(l-f-eMAf)^ « 3{E + eMM), i-e. , 12%, or even 15% 
of the potential events. The probability to miss poten- 
tial events because of the ineffectualness of closed-form 
waveform models will become much more significant in 
searches performed with advanced detectors (see top pan- 
els of Fig.ini). This indicates the need to use more effec- 
tual templates, and to construct template banks having 
a smaller maximal mismatch. 



B. Accuracy. Measurement 

As recalled above, if one wishes to use some waveform 
model not only for detecting the presence of a GW signal, 
but also for measuring the physical information contained 
in the GW signal, one needs this waveform model to sat- 
isfy a stricter condition than the one of effectualness stud- 
ied in the previous section. Namely, the waveform model 
must be such that its inaccurcLcy 1.\hni] 

], Eq. (EZl), 

is smaller than l/p^ff, where poff = p/e denotes what 
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we called the effective SNR (which is larger than the real 
SNR p by a safety factor 1/e > 1), see Eq. Thus, the 
accuracy condition for each model shall be well defined 
after choosing: (i) what is the target signal that one is 
trying to measure; and (ii) a certain level for the effective 
SNR, pos = p/e. As before, we shall choose as fiducial 
target signal the NR-calibrated EOB waveform, EOBnr, 
whereas the effective SNR levels that we shall consider 
are peS = 5, 7, 10, 20. Note beforehand that most of these 
levels are quite lax, and that only the last one is likely 
to be a really useful accuracy requirement. [Indeed, we 
expect typical first observations to have p ^ 10, and, as 
discussed above, adding a safety factor e 1/2 seems to 
be a minimal requirement.] 

Following the same strategy as in the previous Sec- 
tion, we start by studying how the inaccuracy evolves 
with the total mass of the system for a particular mass 
ratio, 4:1, and a particular noise PSD, namely Advanced 
LIGO's. The results are plotted as light gray (red on- 
line) lines in Fig. [5] for the three closed-form models 
we have considered (using EOBnr waveforms as target). 
As expected from our effectualness results (remembering 
the fact mentioned above that the inaccuracy is a non- 
minimized version of - twice - the ineffectualness), we 
find that the PN(f ) model is only accurate for low mass 
systems, while the opposite happens for the phenomeno- 
logical models. For instance, if we assume an effective 
SNR of 10, we obtain that the PN(f) model is only ac- 
curate for describing systems lighter than SMq, whereas 
PhenVl is not accurate for any system we have analyzed 
within M G (3, 500)AfQ and PhenV2 only provides an 
accurate representation for systems heavier than 368Mq. 
In this case, we are finding a very wide inaccuracy inter- 
val, where neither PN(f), nor any of the phenomenologi- 
cal models provide waveforms accurate enough to extract 
faithful information from the GW data. 

We have determined the inaccuracy intervals of the 
three analytical models^° for different mass ratio cases, 
using also different detector noise PSDs and setting dif- 
ferent SNR levels; the results are summarized in Ta- 
ble |TJb). Comparing the right part of this Table ("inaccu- 
racy" ) to its left part ( "ineffectualness" ) we can see how 
stricter the accuracy condition is, compared to the effec- 
tualness one. Typically we are now finding that PN(f) 
waveforms are only accurate up to a few Mq's, leaving 
then a wide gap of at least IOOMq until phenomeno- 
logical models provide an accurate representation of the 
observed waveform. Note in particular that we find in 
nearly all cases that the expected paradigmatic BBHs 
with TOi ^ 7712 ^ 10 — 15Mq cannot be accurately de- 
scribed by any of the closed-form waveforms. In other 
words, we find that, in a large fraction of the parameter 
space, no Fourier-domain closed-form models can be used 



As in Sec. I VI Al their left-most boundaries are determined by 
the PN(f) model and the right-most ones by PhenV2. 



for measurement purposes. 

In order to better visualize the inaccuracy gaps 
that we are finding in parameter space, we repre- 
sent in Fig. |6] (bottom panels) the contour plots cor- 
responding to the inaccuracy levels used in Table HI 
namely 0.04,0.0204,0.01,0.0025, corresponding to pas = 
5, 7, 10, 20. This means that the darkest zone is the re- 
gion where the inaccuracy^^ is larger than 1/5^ = 0.04, 
the next, less dark, zone is the region where I > 1/7'^ = 
0.0204 (but T < 0.04), etc., until the hghtest [but non 
white] zone where I > 1/(20)^ = 0.0025 (but I < 0.01). 
Note how the darkest zones in the bottom figures, where 
the inaccuracy is larger than the very lax level 1/5^ = 
0.04, span (even for initial detectors!) a very wide re- 
gion of parameter space. For advanced detectors, this 
region of "large inaccuracy" can be schematized as: (i) 
if 1 < 9 < 4, aU the systems with lOAfg < M < 75Mq 
belong to this large-inaccuracy zone, while (ii) if <? > 4, 
all systems, whatever be their mass, belong to this large- 
inaccuracy zone^^! 

The main conclusion of this study is that closed-form 
models can be used, for measurement purposes, only 
in two small islands of parameter space: (i) the verti- 
cally elongated white islands on the left of the bottom 
figures, corresponding to small-mass systems (say with 
M < 5Mq), where PN(f) gives an accurate waveform; 
and (ii) the small whitish islands on the right of the bot- 
tom figures, corresponding to some special, large- mass 
systems (with M ^ 200AfQ and q ^ 2), which can be ac- 
curately described by PhenV2(f). In most of the rest of 
parameter space, the inaccuracy of the closed-form mod- 
els is larger than the minimal level required for extracting 
faithful physical information from the GW observations. 

C. Implications of the non- validity regions 

Let us consider in more detail the meaning and impli- 
cations of the presence of extended regions in parameter 
space where closed-form models arc either ineffectual or 
inaccurate (dark, and darkish, regions in Figs. [5] and [7]). 

To have a better grasp of the physics involved in the 
systems behind the abstract (M, v) parameter-space rep- 
resentation, we have overprinted on Figs. [S] and [7] four 
dashed lines which gauge the number of orbital cycles 
(before merger) contributing 99% of the SNR. [The con- 
sideration of a fraction of the SNR of order 99% is nat- 
ural in our discussion where we are ultimately interested 
in effectualnesses of that order or better.] We recall 
that the SNR^ is given by the logarithmic frequency in- 
tegral of the squared ratio r(/) = {hs{f) /hn{f)Y be- 



Note that we plot as inaccuracy domains at some level X the 
region of parameter space where both PN(f) and PhenV2 have 
an inaccuracy larger than I. 

Let us recall that phenomenological models were created using 
NR data only for q < A. 
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tween the (dimensionlcss) effective GW signal ampli- 
tude hs{f) and the (dimensionless) effeetive GW noise 
hnif) = ifSn{f)y^^- The integrand r(/) is maximum 
at some frequency fm and decreases on both sides of /,„. 
For each given detector, each given fraction of SNR (say 
0.99), and each given BBH system, we can then find two 
frequencies, say /i,/2, with /i < /m, < /2, and such 
that/^^^'r(/)dln/= (0.99)^ r{f)d\ii f. Then, we can 
count the number of orbital cycles, say -/Vp.gg, (defined 
by the underlying EOB dynamics) spent by the consid- 
ered BBH system between the GW frequencies /i and 
/2. [When /2 is reached after the merger, one counts 
the orbital cycles between /i and merger.] Actually, the 
choice of the two frequencies is not unique as we can move 
/i and /2 so as to keep the integral jj^ r{f)dlnf fixed. 
But we (numerically) looked for the values of /i and /2 
that minimize the number of orbital cycles spent between 
them. The result is unique, and defines a number of or- 
bital cycles -/Vo.99 which is a function of the parameters 
of the system (i.e. , in our case, of Al and ly) and of the 
choice of detector noise curve. [In most cases, iVg.gg is 
dominated by the number of orbital cycles spent in the in- 
spiral, and is approximately given by the first term (with 
/o fi) of the Newtonian estimate Eq. (|46| . namely 
^^0.99 ^ (7rM/i)~^/3/(647rz/). But this formula does not 
help to know the precise value of iVg.gg because the value 
of the minimizing lower frequency /i is itself a function 
of M, V and the detector PSD, which depends very much 
on the relative location of the peak of the detector sensi- 
tivity curve l/(/i„(/))^, and the merger.] Going from left 
to right, the four dashed lines in Figs. [6] and [7] represent 
the locus of the systems for which A''o.gg(M, ;y) is con- 
stant, and successively takes the values {100, 50, 20, 10}. 
Clearly, for a given z/, iVo.gg(M, i/) will increase as the 
total mass M decreases. This means in particular that 
any system lying on the left of the left-most dashed line 
will have iVg.gg > 100, and similarly for the other lines. 

One should view the level lines A'o.gg ~ const, as 
a (partial) "intrinsic coordinate system" in parameter 
space. That is to say: a level line A'o.gg = const, has 
a more direct meaning, from the data analysis point of 
view, than a coordinate line M = const. These lines al- 
lows one also to make a connection with the capabili- 
ties of NR simulations. Indeed, the current longest NR 
simulations have at most iV ~ 20 orbital cycles before 
merger. Anticipating faster NR simulations, we can hope 
that they will be able to explore in the foreseeable future 
iV ^ 40 orbital cycles, and possibly N ~ 50. On the 
other hand, it will probably be very difficult for NR to 
simulate > 100 orbital cycles. We shall discuss below the 
issue of the minimum number of orbital cycles from NR 
simulations which would be needed to inform us, in the 
most effective way, about the dynamics and radiation of 
BBHs. 

Having in mind the significance of the dashed lines 
in the figures, we can use them to highlight some in- 
teresting physical facts that appear on the inaccuracy 



plots, i.e. the lower plots of Fig. |6l Let us first note 
that these plots exhibit the existence of two "islands of 
accuracy" (i.e. whitish zones), within an "ocean" of in- 
acceptably large inaccuracy (darker zones). [For defi- 
niteness, we shall define this ocean of inaccuracy as the 
union of the two darkest zones.] The (vertically elon- 
gated) island on the left of each plot lies in the small- 
mass region and corresponds to BBH systems for which 
the inspiral PN(f) waveform provides, by itself, an accu- 
rate description of the full EOBnr waveform. The other 
whitish island on the right of each plot (roughly centered 
around ^ 200Afo, g ^ 2) corresponds to BBH systems 
for which the PhenV2(f) waveform accurately agrees with 
the EOBnr waveform. 

• The first fact is that the left islands of accu- 
racy (which we can call the islands of accuracy of 
pure PN waveforms) possess as right boundaries 
{i.e. their "shores" in the dark ocean of inaccu- 
racy) some lines which are approximately parallel 
to the dashed, iVg.gg = const, lines, with large val- 
ues of the constant. For instance, for initial LIGO, 
the right boundary of the left (PN) island is ap- 
proximately located along the iVg.gg = 70, while 
for advanced LIGO it is approximately along the 
A'o.gg = 100 line. For advanced Virgo the value of 
A'o.gg on the boundary is significantly larger than 
100. 

• The second fact, is that the right islands of accu- 
racy (i.e. the islands of accuracy of phenomenolog- 
ical waveforms) do not extend, on their left side, 
beyond A^o.99 = const, lines, with a rather small 
value of the constant, say A'o.gg ~ 5 or so. Notice 
also that the accuracy zone of the phenomenologi- 
cal waveforms seems to be restricted to g < 4. This 
probably reflects the fact that the construction of 
phenomenological models did not correctly incor- 
porate the known PN inspiral phasing, and it is 
consistent with the parameter range of NR wave- 
forms used to build these models [ij. 111] . 

The flrst fact is easily understood. Indeed, the PN(f) 
model does not include any merger signal. Therefore it 
can have a small inaccuracy with respect to the target 
signal (which include the merger) only for systems whose 
SNR is nearly completely given by the inspiral signal. 
The second fact is a priori more surprising, because the 
motivation for the construction of phenomenological sig- 
nals was to approximately describe a waveform compris- 
ing all the needed components of the full signal, namely 
a PN-type inspiral, a NR-described merger, and a lead- 
ing quasi-normal-mode-type ring-down. However, we can 
see in Fig. [5] that the Phen waveforms exhibit significant 
disagreements, both in amplitude and in phase, with the 
target (EOBnr) signal. More quantitatively, if we re- 
member, from the discussion in the Sec. IID above, that 
an inaccuracy smaller than 1/7^ implies that the noise- 
weighted quadratic sum of the fractional amplitude, and 
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phase, difFcrcnccs must be smaller than 1/7 « 0.14, we 
understand upon looking carefully at Fig.[2]that the Phen 
models are too coarse approximations to be 98% accurate 
during most evolutions that will have a significant mix- 
ture of inspiral and merger. The conclusions we can draw 
from these remarks will be discussed in the next section. 

The lack of validity of these analytical models in cer- 
tain regions of the parameter space, will require the use 
of more accurate (and also more time consuming) mod- 
els when performing a search, such as, either full NR 
simulations, or a new and more complex phenomenolog- 
ical model, or EOB waveforms. In particular, the lim- 
ited accuracy of PN(f) that we are observing, specially 
at more extreme mass ratios, will require filling a large 
gap of the BBH evolution between where PN-SPA fails 
and the merger (it may represent several hundreds of or- 
bital cycles, as we have seen in Sec. El), that will probably 
be computationally unachievable for full NR simulations, 
but much more reasonable for semi-analytical methods 
such as the EOB approach. 



VII. CONCLUSIONS 

Our knowledge of the dynamics of coalescing BBH sys- 
tems and their gravitational emission has experienced a 
huge progress during the last 5 — 10 years, both with the 
improvement of analytical methods, namely PN expan- 
sions and EOB approaches, and the breakthroughs that 
occurred in NR. A crucial question at this stage is to 
know whether the recently accumulated PN+EOB+NR 
knowledge is sufficient for the GW data analysis needs of 
the present and upcoming GW detectors. We can clas- 
sify the currently available model waveforms into several 
categories: 

• closed-form, frequency-domain models which are 
only accurate in a limited range of the evolution, 
e.g. PN(f), ring-down models and (NR-fitted) phe- 
nomenological models (which try to extend their 
validity range by combining information coming 
from PN(f), ringdown(f) and some NR simula- 
tions); 

• EOB models, which cover the whole evolution, and 
which are defined in the time-domain by integrat- 
ing ODEs; they can easily incorporate information 
coming either from PN and from NR; 

• full NR simulations, which are necessary for 
describing the nonperturbative physics around 
merger, but can only cover a limited number of 
orbital cycles (< 20) because they are very time 
consuming in terms of computational costs; 

• time-domain PN models (either using closed-form 
expressions or integrating ODEs) , which are limited 
to the inspiral part of the evolution; 



• various types of hybrid models, joining together 
(either in the time-domain, or in the frequency- 
domain) early PN-type waveforms to later NR 
waveforms, so as to cover the whole evolution. 

The time needed to generate a waveform in one of these 
categories varies by many orders of magnitude between, 
say. a closed-form frequency-domain waveform (a frac- 
tion of a second on a single CPU); an EOB waveform 
(several seconds, also, on a single CPU); or, a full NR 
simulation (^ a month on a computer cluster with sev- 
eral hundreds of CPUs). From the data analysis point of 
view, it would be quite useful to be able to generate wave- 
forms as fast as possible [actually, the real issue is to be 
able to access a dense bank of them very fast]. Though 
some trade off between speed and accuracy can be al- 
lowed for, there are minimal accuracy requirements that 
waveforms have to satisfy. This is why our present pa- 
per has focused on a detailed analysis of the accuracy of 
the fastest existing waveforms, namely the closed-form, 
frequency-domain ones. 

In this paper, we have performed an extensive study, 
in the whole parameter space, of the validity of the 
closed-form, frequency-domain models (for non-spinning 
BBHs): i.e. PN(f) and phenomenological models. We 
have studied both their effectualness (detection) and 
their accuracy (measurement) using as target model the 
currently most accurate NR-calibrated EOB waveform. 
We have considered both initial GW detectors and ad- 
vanced ones. 

Our conclusions concerning closed-form, frequency- 
domain waveform models are: 

• When considering initial interferometers, these 
closed-form models are just effectual enough to be 
used for detection. Though their effectualness 

is, indeed, larger than 97% in nearly all the param- 
eter space^'^, it is smaller than 98% in two separate 
islands in parameter space, and smaller than 99% 
in sizable regions of parameter space. See Icf-top 
panel of Fig. [HI and left panel of Fig. [T] In particu- 
lar, notice (on the smaller sub-panels) that most of 
the effectualness is provided by the phenomenolog- 
ical model, except for the very low mass end, and 
also around the crucial M ~ IOM0, q ^ 1 region 
where PhenV2 fails to be 97% effectual. 

• When considering initial interferometers, they are 
not accurate enough to be used for measurement. 
Their inaccuracy, Eq. (|27[) . is, indeed, larger than 
the threshold corresponding to an effective SNR of 
5 (which is an very low threshold) in a huge do- 
main of parameter space. See left-bottom panel of 
Fig. [6] where the darkest zone is the region where 



Apart from a small island near M ~ 20Mq and u = 0.25 in the 
PhcnVl plot, Fig. [7] 
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the inaccuracy is larger than 1/5^. The union of 
the darkest and second darkest zone correspond to 
an inaccuracy larger than 1/7^, which is still a very 
low threshold. 

• When considering advanced interferometers, they 
are neither effectual, nor accurate in large domains 
of parameter space. More precisely, the domains 
of ineffectualness are the extended islands given by 
the dark zones of the two right, upper panels in 
Fig. m (the darkest zone failing to be 97% effectual, 
and the second darkest one failing to be 98% effec- 
tual). Concerning the domains of inaccuracy they 
have become so prominent and extended in the two 
right, lower panels in Fig. [S] (with respect to the do- 
mains of ineffectualness in the plots above them) 
that it is better to visualize them as an ocean of in- 
accuracy containing only two well-separated islands 
of accuracy (namely the whitish zones in the two 
right, bottom panels). In the left island of accu- 
racy, PN(f) provides an accurate model for small- 
mass systems, while, in the right island of accu- 
racy, PhenV2 provides an accurate model for sys- 
tems with total mass M ~ 200Mq, and mass ratio 

Our conclusion is therefore that the current, existing 
closed-form, frequency-domain models fail to be accu- 
rate enough to be used for measurement purposes, even 
with initial GW detectors. On the other hand, if one uses 
jointly the two closed-form models with initial GW detec- 
tors sensitivities, they provide a 97%-effectual coverage 
of the full parameter space, and a 99%-effcctual one in 
an important fraction of parameter space. These effectu- 
alness levels are significantly reduced when considering 
advanced detectors noise PSDs or taking the waveform 
models separately. These results call for the develop- 
ment of more accurate versions of these phenomenolog- 
ical models. In particular, we think that an important 
defect of the VI and V2 versions of the phcnomenological 
models is that they do not correctly incorporate the ana- 
lytically known (PN) phasing behavior during the (early) 
inspiral. One can indeed see on Fig. [2] how this leads to 
large dephasings during the inspiraling phase. In turn, 
these dephasings cause the ineffectualness and inaccuracy 
properties. [See, below, additional results concerning the 
recently constructed V3 version of the phenomenological 
models.] 

In terms of the (partial) intrinsic coordinate system 
defined by the level lines A'o.gg = const., we can summa- 
rize the performance of PhenV2 model saying that it can 
provide an accurate representation of the last A^o.gg ^ 5 



In other words, we can say that the (dark) islands of ineffectu- 
alness of the upper panels have "percolated" among themselves, 
and left only a reverse landscape where the accuracy zones exist 
only as separate islands. 



orbital cycles before merger for g < 3 systems, and an 
effectual representation of signals with no more than 
~ 200 visible orbital cycles within the mass ra- 
tio range 1.5 < q < 8. [PhenVl model is not accu- 
rate enough anywhere and its effectualness is restricted 
to signals with ^ 100 visible orbital cycles in the 

same mass ratio range as PhenV2]. Of course, as we 
consider more sensitive detectors, these upper limits on 
A^o.99 translate into a restriction of a larger portion of the 
parameter space. 

Also, a side result of our study is the determination of 
the maximum frequency, as as function of the mass ratio, 
at which PN(f) waveforms can be 'joined' onto EOBnr 
(in the frequency domain) without undue loss of accu- 
racy. Our results are summarized in Fig.|4]and show wide 
validity ranges for PN(f) up to less than 6 orbits before 
merger for 1:1 and 2:1 systems, but much more worrying 
results when considering q > 4 systems. In particular, we 
find gaps of TV > 20 orbital cycles for peff = 10 [iV > 100 
orbital cycles for p^tf = 20] between the low-end of the 
PN(f) validity range and the merger; which can hardly 
be covered with (computationally intensive) NR simula- 
tions. 

We are fully aware that the analysis that led us to 
all these conclusions depends on our use of the EOBnr 
waveform as a reference model for defining the inaccuracy 
functional Eq. (gT]), Eq. Strictly speaking the dark 

zones in Fig. |6] only represent the regions of parameter 
space where the fractional squared Wiener distance be- 
tween the closed-form models and EOBnr is larger than 
1/7^ « 0.02. However, we think that there is a good 
rationale for considering that our results tell us some- 
thing about the exact inaccuracy of the closed-form mod- 
els, i.e. their fractional squared Wiener distance away 
from the (unknown) exact waveform. Indeed, one should 
keep in mind that the inaccuracy level corresponding to 
the two darkest zones corresponds to fractional ampli- 
tude, and phase inaccuracies (in radians), which are only 
^ 1/7 0.14. By comparison, the NR-calibrated version 
of the FOB formalism that we used here [s^ has been 
proven to agree with some of the current most accurate 
NR simulations with significantly smaller amplitude and 
phase inaccuracies. More precisely, in the q = 1 case, the 
phase disagreement A0 between FOB and the Caltech- 
Cornell data [l^[22| remained smaller than ±0.02 radians 
during the entire 16 orbital cycles of inspiral and plunge, 
which is comparable with the estimated NR errors. At 
the merger and during the ringdown Acj) took somewhat 
larger values (±0.1 radians), but it oscillated around zero 
so that, on average, it stayed very well in phase with 
the NR waveform. As for the fractional amplitude dis- 
agreement |AA|/A it remained smaller than about 0.003 
during the entire inspiral, and only momentarily exhib- 
ited a larger value (0.025) just after the merger. In the 
q = 2 and q = 4 cases, the agreement between EOB 
and the Jena data (published in [3^) was excellent, and 
within the numerical errors during the entire span of the 
simulation (which comprised 17 GW cycles); see Fig. 10 
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in [37|. However, the numerical errors in the latter q = 2 
case were significantly larger than in the q = 1 case, 
but smaller than the ^ 1/7 w 0.14 level used in our in- 
accuracy study (remaining probably smaller than ±0.1 
radians). Let us also mention that an earlier version of 
the EOB formalism has also been succesfully calibrated 
(within NR errors) to accurate NR simulations for sev- 
eral mass ratios q ~ 1,2 and 3 [SSll (however, the 8-orbit 
q = 2, 3 simulations used in |38| cover only part of the 
late inspiral, without the merger). For g > 4, no EOB- 
NR calibrations have been published yet, due to lack of 
accurate and long enough NR waveforms available to us 
in this mass-ratio regime; however, let us note that re- 
cent work [sil HH has accurately confirmed the validity 
of the EOB waveform model in the extreme mass ratio 
limit. 

We have also explored how the results would change 
in the hypothetical case where the low frequency cut- 
off of the advanced detectors would be increased from 
lOHz to 20Hz, since improving the low frequency region 
of the ground-based detectors sensitivity is always the 
most challenging task from the experimental point of 
view. All our results remain almost unaltered, except for 
a tiny improvement in the effectualness of phenomenolog- 
ical models due to a reduction of the number of visible 
orbital cycles. In any case, all the conclusions described 
in this article would remain the same. 

By the time this article was ready for submis- 
sion, a new phcnomenological waveform model for non- 
precessing spinning BBH systems was published by San- 



tamaria et al. l93|. This model (that we shall call 



"PhenV3") represents an improvement over the previ- 
ous PhenV2 model. This new PhenV3 model aims at 
describing the whole BBH coalescing process by join- 
ing a state-of-the-art PN description of the early inspi- 
ral (which was lacking in VI and V2) to separate phe- 
nomcnological representations of a "pre-merger signal" , 
and a "ring-down" one. In view of the improved na- 
ture of this PhenV3 model, we have decided to complete 
our analysis by studying its effectualness and accuracy 
performances (i.e. the analogs of Figs. [6] and [7] above). 
Our results for the PhenV3 waveform^'^ arc displayed in 
Fig. [51 As this new model now incorporates a 3.5 PN- 
accuratc description of the early inspiral, we directly plot 
in Fig. [8] the incffectualness and inaccuracy (with respect 
to EOBnr) of PhenV3(f) alone, without computing the 
infimum with respect to PN(f). If we focus our attention 
on the g < 4 region (where the phcnomenological param- 
eters were fitted), we observe that the PhenV3 model is 
> 98% effectual for initial detectors and > 97% effectual 
when considering advanced detectors noise PSDs. The 
effectualness is higher for massive systems, M > 50Mq, 



The authors of [o^ l kindly let us know about the existence of 
a typo in their Eq. (5.12): the 'a' parameter of the Lorentzian 
function should read Sifsn/Q instead of 52Q- Here, we have used 
the corrected expression when generating our results. 



(where only the last few cycles of the evolution are ob- 
servable) and reaches minimum values in the [astrophys- 
ically relevant] region around M ~ IOMq (where there 
are still many observable cycles but PN theory has al- 
ready started to fail). Notice that these results are con- 
sistent with the overlaps plotted in Fig. 10 of [l^l, despite 
the fact that here we do not maximize the effectualness 
over any spin parameter, as we have restricted ourselves 
to non-spinning systems. The fact that the ineffectual- 
ness (as well as the inaccuracy) is maximum around q ~ 2 
and M ^ 20A1q, is compatible with the fact that, when 
plotting the q = 2 analog of Fig. [U for the PhenV3/E0B 
comparison, we observed the existence of significant de- 
phasings, of order ^ 0.5 rad, around fM ~ 8x 10~^. Sim- 
ilarly to what was observed with PhenVl and PhenV2 
models, the PhenV3 effectualness becomes much worse 
for systems with q > 4. In addition, it is important 
to note that the accuracy of the PhenV3 model is only 
slightly improved with respect to the PhcnV2 one. Actu- 
ally, Fig. [5] shows that the PhenVS model cannot be used 
for measurement purposes (accuracy) in most of the pa- 
rameter space. This shows that the conclusions drawn 
above from a study of the VI and V2 models hold also 
for the PhenV3 model. 

Also, while writing up our results for publication, an- 
other paper related to our work appeared on the archives. 
Indeed, the article |94| . presented another study similar 
in spirit to the one of our Sec. V. However, the con- 
text, and details of their study are quite different from 
ours: they are using as target waveform a (spinning) 
PNU/^j^^jNR hybrid waveform (where the starting fre- 
quency /start of their NR waveform covers 6-10 orbital 
cycles before merger) and they compute the faithfulness 
between this target and a hybrid of the form PNU/qNR, 
where the junction frequency /o is larger than /start- In 
other words, they are monitoring the loss of effectualness 
in PNUNR hybrids as one includes less and less NR or- 
bital cycles in the hybrid (starting from 6-10 NR orbital 
cycles). 



VIII. FUTURE WORK 

In view of our results, we think that, especially for 
future searches in advanced detectors data, more effec- 
tual and more accurate model waveforms will need to be 
considered in order not to miss a significant fraction of 
the hard-won potential detections and, to extract reli- 
able physical data out of the observations. Though some 
improved (more effectual) phcnomenological waveforms 
might be useful for detection purposes, we think that the 
NR-calibratcd EOB waveforms are our current "best bet" 
models for building waveforms that are both effectual and 
accurate. More effort should be put into improving the 
EOB waveforms as well as into improving their imple- 
mentation in data-analysis platforms. Here is a list of 
suggestions towards this aim. 

• More work is needed to refine the results of Sec. V, 
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FIG. 8: Same plots as in Fig. (6] but considering the recently published PhenV3 [l^l model. Since it already includes all the PN 
contributions in the low frequency region, here we consider the phenomenological model alone. We recall that the top panels 
plot the regions of the parameter space where the ineffectualness is larger than {3%, 2%, 1%, 0.45%} going from the darkest to 
the lightest areas; whereas the bottom panels show the regions where the inaccuracy is laxger than 1/Pcff i where the effective 
SNR is pcff = {5, 7, 10, 20}, when going from the darkest to the lightest region. 



displayed in our Fig. 2] concerning the maximum 
junction frequency M fo at which one can smoothly 
join (in the frequency-domain) an early PN(f) 
waveform to a later EOB(f) one without losing 
more than some given accuracy. Our results sug- 
gest that it is ill-advised to continue the common 
practice of defining (either in the time-domain, 
or the frequency-domain) hybrid PNU/(,NR wave- 
forms by joining an early PN waveform to a later 
NR one, at a frequency fo near the starting fre- 
quency of the NR waveform. Such a procedure is 
probably justified (without losing too much accu- 
racy) only for NR waveforms containing more or- 
bital cycles before merger than the numbers dis- 
played in the last columns of the Table contained 
in Fig. 21 This would make it practically impossi- 
ble to build accurate PNUNR waveforms for mass 
ratios g > 4. 



better quality of the resummed expressions^^ used 
in the BOB formalism will allow one (if one insists 
on doing so) to join an early EOB waveform to a 
later NR waveform around a frequency /o corre- 
sponding to significantly less orbital cycles before 
merger than the numbers displayed in the Table 
contained in Fig. |4l We leave to a future work [7l| 
a discussion of the minimum number of NR orbital 
cycles needed to improve the NR-calibation of EOB 
to the level of accuracy needed for the analysis of 
the future GW data. 

• Independently of the quantitative issue of this min- 
imum number of NR orbital cycles, we think that 
it is more efficient to use accurate NR data to di- 
rectly inform the EOB formalism by calibrating 
some unknown higher-order theoretical parameters, 
and then to use the resulting NR-calibrated EOB 
waveforms as best-bet models (rather than passing 



• We think, however, that most of the accuracy loss 
behind Fig. 2] is due to the poor quality of (any 
of) the non-resummed PN expansions^ and that the 



The effectiveness of the resummation methods used in the most 
recent version of EOB is exemplified, for instance, in Fig. 1 of 
[ToII for the u -t^ I case, and in Fig. 4 of [s^l for the u = 1/A case. 



28 



through the intermediate step of building hybrid 
EOBUNR waveforms). These EOBnr waveforms 
could then be used (if wished) as basic material 
for defining more accurate closed-form frequency- 
domain phenomenological-type waveforms. 

• However, we think that it should be possible to 
set up a (multi-point) interpolation technique us- 
ing a relatively sparse bank of (pre-computed) 
EOBNR(f) waveforms to compute, in a numeri- 
cally fast way, any wished EOBNR(f) waveform, 
with the accuracy required for measurement pur- 
poses. A detailed analysis is needed to estimate 
the CPU time needed, at some required inaccuracy 
level I, by such an interpolation technique, which 
would then obviate the need for having at hand 
a closed-form, frequency-domain model. [However, 
improved closed-form models might still be useful 
for detection purposes.] 

Let us note finally that several ongoing works might 
soon allow one to define more accurate analytical de- 
scriptions of BBH systems. First, a recent joint effort 
between several numerical relativity and analytical rel- 
ativity groups [zll has begun to make use of TeraGrid 
resources in order to produce longer and more accurate 
simulations of comparable-mass BBH systems. Second, 
the case of extreme mass-ratio BBH mergers is amenable 
to separate Regge- Wheeler- Zerilli- type studies which in- 



form us about the structure of the waveform in the z/ ^ 1 
limit [H, iH, illil- Third^rogress in Gravitational 
Self Force (GSF) studies [8l,|83| has recently allowed one 
to get a new, accurate source of information about the 
strong-field behavior of some of the basic ingredients of 
the EOB formalism [H, [s^ . We note also that our cur- 
rent study has only considered the leading quadrupolar 
waveform {i = 2,m — ±2), and non-spinning systems. 
The EOB formalism provides (resummed) waveforms for 
all multipolarities m [t^I so that it will be straight- 
forward to extend our analysis to the full waveform (at 
least for the comparison with the corresponding higher- 
multipolarity FN waveform). The EOB formalism has 
been defined for spinning BBH systems [H EO, 
so that our study can also be extended to the case of 
spinning BBH waveforms. 
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